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SECTION  I 


INTRODUCTION 

The  Department  of  Transportation  Transportation  Systems  Center  (DOT  TSC) 
is  developing  a Multi-Modal  Transportation  System  Simulation  with  real-time 
capability  for  simulating  the  command  and  control  functions  related  to  various 
transportation  systems  (Ref.  1).  The  first  phase  of  this  effort  is  producing  an  Ad- 
vanced Air  Traffic  Control  (ATC)  Simulation  Laboratory  to  define  and  evaluate 
the  most  effective  role  of  the  controller  in  advanced  generation  ATC  systems,  and 
to  assess  quantitatively  various  control  strategies  and  automated  techniques  in  these 
systems.  The  objective  is  to  provide  and  exercise  a system  of  interactive  experi- 
ments to  verify  concepts  of  advanced  ATC  systems  using  fast-time  and  real-time 
computer  simulations. 

The  types  of  experiments  to  be  investigated  with  the  TSC  simulation  include, 
for  example:  strategic  control  concepts  where  aircraft  are  assigned  four-dimensional 
(4-D)  route-time  profiles  developed  on  the  ground  prior  to  flight  (Refs.  2 - 4); 
sequencing,  spacing  and  flow  control  concepts  for  regulating  the  traffic  volume 
in  the  congested  terminal  area  (Refs.  5 - 8);  and  separation  assurance  concepts 
such  as  intermittent  positive  control  (IRC)  and  collision  avoidance  systems  (Refs. 

5,  9,'  10). 

The  ATC-oriented  program  is  planned  to  extend  over  a period  of  several 
years.  Since  continuing  investigations  are  expected  to  provide  new  ATC  knowl- 
edge, the  programs  created  for  the  simulation  are  being  made  extremely  flexible 
and  easy  to  expand  and  change.  In  addition,  as  new  experiments  are  defined, 
the  simulation  programs  must  be  adapted  to  provide  the  user  with  the  necessary  tools 
for  his  purpose. 
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The  present  report  presents  the  results  of  a study  conducted  by  Aerospace 
Systems,  Inc.  (ASI)  to  provide  more  realistic  error  models  for  the  TSC  simulation. 

The  original  A $ 1 objective  was  to  evaluate  selected  control  strategies  for  satisfying 
4-D  guidance  requirements  and  to  determine  the  effects  of  the  more  important 
system  errors  by  simulation.  However,  it  became  apparent  that  the  error  models 
themselves  could  easily  dominate  the  results  of  such  an  analysis.  Consequently, 
the  ASI  effort  concentrated  on  developing  more  sophisticated  models  of  the  dominant 
functional  error  sources  to  support  the  ATC  systems  evaluation  experiments  planned 
for  the  TSC  simulation. 

The  next  section  of  the  report  presents  a broad  overview  and  background 
of  the  present  ATC  simulation.  The  requirements  for  certain  improvements  in  the 
simulation  are  mentioned  and  the  general  philosophy  followed  in  designing  the 
revised  models  is  discussed.  The  next  four  sections  present  the  revised  error  models 
for  the  simulation:  Section  3 describes  the  motions  of  the  aircraft;  Section  4 presents 
the  Air  Data  System  errors;  Section  5 describes  the  Navigation  System  Models; 
and  Section  6 presents  the  Ground  Surveillance  and  Radar  Tracking  models.  The 
nomenclature  used  in  these  models  is  defined  in  the  front  of  the  report.  This  nomen- 
clature and  the  corresponding  units  are  consistent  with  that  used  in  the  analysis 
and  operations  of  ATC  systems.  In  Section  7,  various  features  of  the  revised  models 
are  discussed,  and  some  comparisons  are  made  with  the  previous  models.  The 
significant  conclusions  of  the  study  and  suggested  applications  for  the  revised 
models  are  presented  in  Section  8.  Finally,  a comprehensive  list  of  References  , 
and  a Bibliography  appear,  upon  which  much  of  the  modeling  was  based. 
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SECTION  2 


BACKGROUND 

This  section  presents  some  background  for  the  development  of  the  ATC 
error  models  described  in  the  subsequent  sections  of  the  report.  The  present  TSC 
ATC  simulation  is  reviewed  and  the  need  for  certain  revisions  is  indicated.  The 
philosophy  adopted  for  developing  the  revised  models  is  discussed  briefly. 

2.1  ATC  SYSTEM  MODEL 

A general  block  diagram  of  an  air  traffic  control  system  as  it  applies 
between  a single  aircraft  and  controller  is  shown  schematically  in  Figure  1 . 

The  target  aircraft  submits  a flight  plan  request,  F , to  the  sector  controller  who 
in  turn  issues  a set  of  commands  in  the  form  of  a clearance,  Fc,  to  the  aircraft 

•t. 

for  execution.  The  clearance  would  normally  consist  of  a route  for  the  aircraft 
to  follow  and  an  altitude  and  airspeed  to  maintain  to  prevent  conflict  with  any 
other  controlled  aircraft.  Under  vectoring  control  in  the  terminal  areas  the  clearance 
would  specify  a heading  to  maintain  instead  of  the  route.  For  advanced  control 
concepts,  the  clearance  might  specify  a groundspeed  or  the  times  at  which  to 
cross  certain  fixes.  For  collision  avoidance,  the  clearance  might  call  for  climb, 
dive  or  turn  maneuvers,  followed  by  a command  to  resume  normal  navigation. 

The  air  data  system  provides  the  aircraft  commander  with  information 
about  the  heading,  altitude,  altitude  rate,  airspeed  and  Mach  number  of  his  air- 
craft. The  navigation  system  provides  information  about  the  location  of  the  air- 
craft relative  to  the  desired  ground  track.  Both  the  air  data  and  navigation  systems 
information  contain  some  errors.  The  aircraft  commander  uses  the  information  pro- 
vided by  these  systems  to  control  the  velocity  vector  V , of  the  aircraft  relative 
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STRATEGY  NOT  KNOWN 
ERRORS  NOT  MODELLED 


STRATEGY  KNOWN 
ERRORS  MODELLED 


Figure  1 . General  Air  Traffic  Control  System  for  a Single  Aircraft. 
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to  the  air  mass,  in  a manner  which  he  believes  will  result  in  adherence  to  his 
clearance.  In  so  doing,  he  uses  control  strategies  which,  fortunately,  are  fairly 
standard  and  repeatable. 

The  groundspeed  of  the  aircraft  V , is  the  vector  sum  of  the  airspeed  and 

9 

the  wind,  V . integrates  into  aircraft  position  relative  to  the  ground,  P.  The 
positions  of  target  aircraft  are  sensed  by  the  surveillance  system  and  the  information 
is  made  available  to  the  sector  controller.  The  sector  controller  also  has  some  knowl- 
edge  of  the  wind,  , through  the  weather  observations  and  forecasts.  The  sur- 
veillance information  is  used  to  effect  changes  in  the  clearances  given  to  the  target 
aircraft. 

The  net  result  of  various  advanced  control  strategies  which  might  be  used 
by  the  sector  controller  is  now  not  known,  nor  is  the  effect  of  system  errors  on 
the  operation  under  any  given  strategy.  The  overall  goal  of  the  TSC  effort  is  to 
be  able  to  establish  these  relationships  by  simulation  of  the  system. 

2.2  ATC  SIMULATION 

The  ATC  simulation  laboratory  consists  of  a PDP-10  digital  computer  facility 
and  a display  system,  which  uses  a DDP-516  computer  coupled  to  an  ADDS/900 
display  unit  (Ref.  1).  The  simulation  uses  the  PDP-10  as  the  ground-based  data 
acquisition  and  command  and  control  system,  with  the  DDP-516  computer  as  the 
display  processor,  (Fig.  2).  Basically,  the  functions  that  are  performed  by  the 
PDP-10  computer  are:  (1)  simulation  of  aircraft  traffic  including  models  for  the 
pilot,  aircraft,  terminal  area,  and  environment;  (2)  control  of  traffic  with  automation 
and  decision-making  algorithms;  (3)  manipulation  of  data  files  and  control  to 
update  the  display  list;  (4)  scenario  generation  for  each  simulation;  (5)  data  re- 
cording and  reduction;  (6)  control ler/pi lot-aircraft  command  processing;  and 
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Figure  2.  Air  Traffic  Control  Simulation. 

(7)  flight  strip  printing.  The  functions  performed  by  the  display  processor  are: 

(1)  generation  of  synthetic  background  video  maps;  (2)  display  file  maintenance; 

(3)  display  refreshing:  (4)  I/O  handling  for  data  entry  devices;  and  (5)  controller 
command  execution. 

The  simulation  software  has  four  basic  elements  (Ref.  1):  traffic  generation, 
sky  and  terminal  environment,  pilot  module,  and  control  module.  The  traffic 
generator  module  contains  the  driving  mechanism  of  the  model  and  introduces  flight 
characteristics,  update  timing  and  sequencing  characteristics.  A scenario  generates 
aircraft  arrivals  and  departures  and  introduces  all  of  the  circumstances  required 
to  simulate  special  events  such  as  airborne  equipment  failures,  runway  closures 
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and  reversals,  and  missed  approaches.  The  environment  module  includes  flight 
characteristic  limits,  control  response,  aircraft  errors,  and  a geographical  section 
that  contains  position  and  route  definitions.  The  pilot  module  accepts  pilot  commands 
for  heading,  speed,  and  altitude  and  then  implements  these  commands  to  the  air- 
craft. The  intent  of  the  pilot  module  is  not  to  model  a pilot  but  to  allow  an  in- 
dividual to  interface  the  audio  commands  of  the  controller  to  the  digital  input  re- 
quired by  the  PDP-10.  The  controller  is  provided  with  a realistic  display  of  traffic 
and  can  issue  commands  to  a pilot  through  an  audio  link  or  a simulated  data  link. 
Details  of  the  simulation  program.  Route  Oriented  Simulation  System  (ROSS), 
are  contained  in  References  1 1 and  12.  A brief  review  of  the  error  models  in  ROSS 
is  presented  next. 

The  wind  environment  was  represented  by  layers  of  constant  winds.  The 
wind  speed  and  direction  were  constant  within  each  layer,  and  changed  instanta- 
neously between  two  layers.  No  vertical  winds  were  included.  At  ASI's  suggestion, 
ROSS  was  modified  to  permit  wind  speed  and  direction  to  be  continuous  linear 
functions  of  altitude. 

The  aircraft  altitude,  speed  and  heading  are  updated,  if  necessary,  as 
linear  functions  of  time.  The  rate  of  climb,  rate  of  descent,  acceleration,  decel- 
eration and  rate  of  turn  are  all  constant  parameters  for  a particular  aircraft.  Eaqh 
parameter  is  chosen  from  a normal  distribution  with  a specified  mean  and  standard 
deviation  appropriate  for  the  aircraft  type.  No  dynamic  response  is  included  in 
the  motion  module;  i.e.  there  is  no  lag  in  responding  to  a command,  and  no  over- 
shoots of  the  commanded  value  are  allowed.  For  example,  if  the  altitude  change 
overshoots  the  desired  altitude,  the  exact  altitude  is  set  equal  to  the  desired  altitude. 
However,  pure  time  delays  can  be  introduced  to  represent  communications  delays, 
for  example.  Fixed  instrument  bias  errors  can  be  represented  by  selecting  the  com- 
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manded  altitude,  heading  or  airspeed  from  specified  random  distributions. 

The  new  airspeed  and  wind  components  are  added  to  obtain  groundspeed  com- 
ponents. The  aircraft  position  is  then  updated  by  extrapolation  from  its  previous  value, 
using  the  new  groundspeed  components.  If  the  aircraft  is  being  observed  by  radar, 
the  range  and  bearing  errors  (each  defined  by  a mean  and  standard  deviation)  are 
determined,  and  the  noisy  position  coordinates  relative  to  the  radar  site  are  calculated. 
When  a navigation  station  is  being  used,  a first-order  noise  process  (with  a 0.03  sec 
time  constant)  produces  sequentially-correlated  range  and  bearing  errors,  which  are 
converted  to  noisy  position  estimates. 

A homing  mode  of  navigation  is  provided  in  which  the  aircraft  flies  directly 
to  a specified  position.  If  the  aircraft  is  more  than  two  miles  from  the  homing  position, 
the  heading  to  that  position  is  computed.  To  account  for  crosswinds,  the  drift 
angle  which  would  result  if  the  aircraft  flew  on  the  heading  to  the  homing  position 
is  calculated.  This  drift  angle  is  then  applied  as  a correction  to  obtain  a desired 
heading  compensated  for  the  wind  effects.  When  the  aircraft  is  less  than  two  miles 
from  the  homer,  a fixed  heading  Is  commanded. 

2.3  IMPROVED  ERROR  MODELS 

The  ROSS  program  was  developed  by  using  relatively  simple  kinematics 
to  model  the  target  aircraft  and  by  neglecting  certain  of  the  errors  in  the  navigation, 
air  data  and  surveillance  systems.  In  general  there  is  a trade-off  to  be  made  between 
the  complexity  of  the  model  representing  the  system  and  the  cost  in  computer  pro- 
gramming and  execution  time.  As  one  needs  to  know  better  the  high  frequency 
response  of  the  system,  it  is  necessary  to  model  the  system  dynamics  in  greater 
detail.  At  the  present  time  the  ROSS  program  generally  provides  unrealistic  dynamic 
response  characteristics  at  the  higher  frequencies,  i.e.  those  frequencies  having 
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periods  below  30  seconds.  The  present  simulation  permits  instantaneous  changes 
in  attitude  and  heading  rate,  for  example.  This  would  require  a much  faster  aircraft- 
pilot  response  than  is  physically  possible.  Consequently,  the  dynamic  response 
characteristics  need  to  be  modified  in  the  high  frequency  spectrum  or,  equivalently, 
for  response  times  below  about  30  seconds. 

There  are  several  reasons  for  wanting  to  extend  the  capability  of  ROSS  to 
represent  high  frequency  effects  more  accurately.  One  of  these  is  to  evaluate  various 
collision  avoidance  systems,  including  intermittent  positive  control.  Typical  warning 
times  for  such  systems  are  less  than  30  seconds  and  separation  distances  will  be  under 
a few  thousand  feet.  Another  application  of  interest  is  the  4-D  guidance  system, 
using  clearances  which  require  the  aircraft  to  hit  navigation  fixes  at  specified  times. 
The  dynamic  effect  of  the  wind  and  the  dynamic  errors  in  the  navigation  system  are 
expected  to  have  a strong  influence  on  how  this  control  strategy  can  be  used.  Another 
interesting  application  is  to  investigate  strategies  for  sequencing  and  spacing.  It  has 
been  estimated  that  timing  errors  of  better  than  4 seconds  are  necessary  to  achieve 
maximum  runway  utilization  (Ref.  10).  The  simulation  would  have  to  include 
all  effects  of  that  order  of  magnitude  in  order  to  assess  the  various  control  strategies. 

The  purpose  of  this  work  has  been  to  expand  the  models  to  the  right  of 
the  dashed  line  in  Figure  1 to  include  important  effects  having  response  times 
longer  than  one  or  two  seconds.  The  philosophy  has  been  to  introduce  accuracy 
and  realism  to  the  greatest  extent  practicable  while  minimizing  the  requirement 
for  extensive  computation.  In  most  cases  the  target  dynamics  depend  strongly 
on  the  control  strategy  used  by  the  aircraft  commander  in  following  his  clearance. 
Fortunately,  pilot  procedures  for  instrument  flight  are  well  standardized  by  FAA 
guidelines  and  regulations.  Every  effort  has  been  made  to  have  the  analytic  models 
reflect  as  nearly  as  practicable  the  actions  taken  by  aircraft  commanders  in  the 
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real  world.  The  models  have  been  developed  by  pilots  familiar  with  instrument 
operations,  using  their  judgement  where  necessary  to  estimate  the  magnitude  of 
system  parameters.  It  has  been  assumed  that  all  clearances  are  executed  upon 
receipt;  communications  delays  in  receiving  clearances  have  not  been  modeled, 
but  could  be  included  by  using  the  pure  delays  presently  available  in  ROSS. 

The  dynamics  of  the  aircraft  have  been  modeled  on  the  basis  of  physical 
laws  to  the  extent  practicable.  Dynamical  modes  of  the  aircraft  which  would  be 
suppressed  by  automatic  or  human  pilot  damping  have  been  neglected.  Errors 
in  the  air  data  system  are  modeled  as  instrument  biases.  The  navigation  system 
errors  are  modeled  separately  for  the  kind  of  navigation  system  under  consideration. 
Errors  and  loss  of  information  due  to  both  geometry  and  propagation  are  included. 
Where  necessary,  errors  are  time-correlated  in  order  to  provide  realistic  system 
response.  The  surveillance  system  includes  the  same  type  of  sensing  errors  in  addition 
to  the  general  algorithms  used  in  the  third-generation  ATC  system  tracking  software. 

A note  of  caution  is  appropriate  regarding  the  units  used  in  the  error 
models.  The  primary  input  and  output  quantities  (e.g.  range,  altitude,  airspeed, 
rate  of  climb)  have  been  expressed  in  the  standard  terms  used  in  the  ATC  environ- 
ment (e.g.  nautical  miles,  feet,  knots,  feet  per  minute).  Intermediate  variables, 
however,  use  feet  and  seconds  as  the  basic  length  and  time  units.  Consequently, 
the  reader  must  be  careful  to  avoid  confusing  these  terms.  In  several  places  the 
models  contain  specific  numerical  values  to  convert  quantities  having  dissimilar 
units . 
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SECTION  3 


TARGET  DYNAMICS  MODEL 

The  target  dynamics  mode!  provides  the  means  of  representing  the  motion 
of  each  aircraft  (target)  in  the  simulated  ATC  system.  The  development  of  appro- 
priate target  dynamics  models  involves  a tradeoff  between  realism  and  computational 
requirements.  If  unlimited  simulation  resources  were  available,  the  ideal  situation 
could  be  a complete  six-degree-of-freedom  simulation  of  each  aircraft  in  the  system, 
including  realistic  aerodynamics,  engine  characteristics  and  a full  cockpit. 
Obviously,  this  is  not  a feasible  alternative,  and  would  not  be  required  for  the 
experiments  being  proposed  for  the  TSC  simulation. 

On  the  other  hand,  the  motion  update  model  presently  existing  in  the 
simulation  does  not  include  any  dynamic  response  of  the  targets;  e.g.  a commanded 
change  in  altitude  is  immediately  followed  by  a constant  rate  of  climb  until  the 
new  altitude  is  precisely  attained  (perhaps  with  some  steady-state  bias  error  from 
the  altimeter).  The  target  dynamics  model  presented  in  this  section  will  provide 
the  simulation  user  with  a compromise  between  realism  and  computational  simplicity. 
The  target  dynamics  model  has  also  been  developed  to  be  as  general  as  possible, 
so  as  to  require  minimum  changes  to  represent  the  entire  spectrum  of  aircraft  using 
the  ATC  system.  The  representation  of  a single  engine  general  aviation  aircraft 
or  a military  fighter  is  determined  by  the  selection  of  an  appropriate  set  of  per- 
formance parameters  and  response  time  constants. 

Figure  3 shows  a general  block  diagram  of  the  target  dynamics  model . 

The  inputs  include  commanded  altitude  (h  ),  commanded  airspeed  (V  ) or  Mach 

C b 

number  (M  ),  and  commanded  heading  (\|r  ) or  track.  Disturbances  entering  the 
c c 

model  include  winds,  air  data  system  errors  and  navigation  system  errors.  The  outputs 
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AIRCRAFT  PARAMETERS 


WINDS,  AIR  DATA  8 
NAVIGATION  SYSTEMS 
OUTPUTS 


Figure  3.  General  Block  Diagram  of  Target  Dynamics  Model. 


are  the  target's  true  altitude  (h),  airspeed  (V^),  Mach  number  (M),  heading  (1)1), 

ground  speed  (V  ),  ground  track  angle  (x),  and  horizontal  position  coordinates 
9 

(x,y).  Individual  elements  of  the  target  dynamics  model  are  discussed  in  detail 
below  according  to  the  commanded  input.  Nomenclature  is  defined  at  the  beginning 
of  the  report . 


3.1  ALTITUDE  COMMAND 

A block  diagram  of  the  altitude  command  loop  is  shown  in  Figure  4. 

The  box  labelled  "cl imb/descend  logic"  generates  an  altitude  rate  command  based 
on  the  current  error  e^  between  the  commanded  and  indicated  altitudes.  If  this 
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13 


Figure  4.  Block  Diagram  of  Altitude  Command  Loop. 


error  is  greater  than  1,000  ft,  the  best  rate  of  climb  or  descent  (HDP  or  HDN) 
is  commanded  (Ref.  13).  The  parameters  HDP  and  HDN  are  functions  of  the  air- 
craft type  and  the  true  altitude.  When  the  aircraft  is  within  1,000  ft  of  the  com- 
manded altitude,  the  commanded  rate  is  500  ft/min.  The  altitude  tolerance  is 
shown  as  50  ft;  i.e.  when  the  altitude  error  is  less  than  50  ft,  no  corrective  action 
is  taken . 

The  commanded  altitude  rate  h^  is  compared  with  the  indicated  h|,  and 
any  difference  is  used  to  generate  a suitable  vertical  acceleration.  The  altitude 
rate  feedback  loop  represents  a first-order  lag  response  with  a time  constant 
t^,  and  having  constraints  (GMXP  and  GMXN)  on  the  magnitude  of  the  vertical 
g-loading.  The  output  of  the  box  labelled  "G  limiter"  is  the  true  vertical  accelera 
tion,  which  is  integrated  once  to  give  true  altitude  rate,  and  again  to  give  true 
altitude.  Instrument  errors  are  introduced  in  the  air-data  system  which  provides 
the  indicated  values  hj  and  hj  . 

With  appropriate  values  of  the  parameters  HDP,  HDN,  GMXP,  GMXN, 
and  t^,  the  model  provides  a reasonably  accurate  dynamic  representation  of  the 
target  response  to  an  altitude  command.  The  values  of  the  parameters  given  later 
in  this  section  (Table  1)  are  for  normal  enroute  or  terminal  operation.  They  do 
not  adequately  represent  situations  such  as  collision  avoidance  or  landing  approach, 
when  response  times  would  be  shorter  and  when  normal  constraints  might  be  ignored. 

The  altitude  command  model  depicted  in  Figure  4 is  based  on  normal 
control  by  a typical  pilot.  However,  no  two  pilots  and  no  two  aircraft  will  ever 
have  identical  characteristics.  For  added  realism,  the  various  model  parameters 
should  be  selected  from  an  appropriate  random  distribution  provided  by  the  user. 

In  this  way,  the  behavior  of  each  target  in  the  simulation  will  be  unique.  In 
certain  situations,  particularly  landing  approach,  the  cl imb/descent  logic  might 
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require  special  modification.  For  example,  the  rate  of  descent  during  final  approach 
should  be  a function  of  groundspeed  in  order  that  the  Minimum  Descent  Altitude 
(MDA)  or  Decision  Height  (DH)  be  reached  at  the  correct  distance  from  the  runway. 

The  equations  of  motion  of  the  model  can  be  integrated  analytically  for 
closed-form  solutions,  since  h^  will  be  piecewise  constant.  If  the  "g  limits" 
are  not  in  effect,  the  altitude  rate  and  altitude  are  given  by: 

h(t)  = J_  <|hc(ti)  - e^(t-|)  exp[-  (t  - t1)K^/xh]l  (1) 

Kh 

h(t)  = h(t1)  + _L-  |hc(t1)[t  - t,]  - Th[fi(t)  - h(t  )]}  (2) 

60  K ^ ' 

where  t^  is  the  time  of  the  last  update  and  is  a scale  factor  error  in  the  air 
data  system  (see  Section  4.2).  When  the  "g  limit"  is  in  effect,  the  equations 
are  even  simpler: 

h(t)  = h(tj)  + 60  h(tj)[t  - tj]  (3) 

h(t)  = h(t i ) + h(t -j )[ t - t^/60  + h(t -j )[t  - 1 1 ] 2/ 2 (4) 

The  latter  equations  could  be  used  when  the  g limit  is  not  in  effect, 
thereby  introducing  a maximum  error  of  less  than  100  ft/min  in  li(t).  The  step 
size  used  in  the  simulation  (t  - t-j)  is  a constant  which  is  generally  set  at  1 second 
in  the  present  version.  The  exponential  function  in  Equation  (1)  could  be  approx- 
imated by  the  first  three  terms  in  the  series  expansion  to  conserve  computer  time. 

This  simplification  would  result  in  an  error  of  about  one  percent  or  less.  If  the  first 
two  terms  of  the  series  are  used.  Equation  (3)  Is  obtained. 

3.2  AIRSPEED  COMMAND 

In  response  to  an  airspeed  command,  the  pilot  will  adjust  the  thrust 
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and/or  drag  to  provide  the  desired  longitudinal  acceleration  or  deceleration. 

However,  the  resulting  acceleration/deceleration  acts  upon  the  ground-referenced 
(inertial)  velocity  of  the  aircraft,  which  is  seldom  the  same  as  the  air-referenced 
velocity  due  to  the  effects  of  winds.  Similarly,  lateral  accelerations  in  turning 
maneuvers  will  also  affect  the  airspeed  response  as  a result  of  crosswinds.  This 
coupled  airspeed-heading  response  will  be  discussed  next. 

Figure  5 depicts  the  situation  of  an  aircraft  flying  through  a wind.  The  wind 
magnitude  is  V^,  blowing  from  the  magnetic  direction  9^.  The  aircraft  heading  in  the 


Y 


X (EAST) 


Figure  5.  Geometry  and  Nomenclature  for  an  Aircraft  Flying 
through  a Wind . 
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air  mass  is  ^ , and  the  airspeed  is  V^.  The  resulting  groundspeed  is  V and  the 

ground  track  angle  is  X-  The  longitudinal  acceleration  A applied  to  the  aircraft 

x r 

is  proportional  to  the  difference  between  the  thrust  and  the  drag,  both  of  which  act 
parallel  to  V in  coordinated  flight.  The  lateral  acceleration  A results  from  a co- 

a y 

ordinated  bank  which  rotates  the  aerodynamic  lift  vector  about  V . The  equations 
of  motion  in  ground-referenced  coordinates  are: 


V - A cos  (if  - x)  - A sin  (il i - x) 
g x y T 


V X - A sin  (if  = X)  + A cos  (if  - x) 

g x y ' * • 


By  geometry,  the  airspeed  and  heading  are: 


V = [ V 2 + V 2 + 2 V V cos  (9  - X)] 1/2 

a g w 9 w v w ' 


if  - X + sin  [(Vw/Vq)  sin  (0w  - X)] 


(5) 

(6) 

(7) 

(8) 


Note  that  if  there  is  no  wind  (V  = 0),  V = V and  if  = X. 

w a g 

Figure  6 presents  the  general  block  diagram  for  the  coupled  aircraft 
response  to  air  speed  and  heading  commands.  The  heading  command  loop  shown 
in  the  lower  part  of  the  diagram  will  be  discussed  later  in  this  section. 

In  the  airspeed  command  loop,  an  error  between  the  commanded  and 
indicated  airspeed  ey  in  excess  of  2.5  knots  results  in  a longitudinal  acceleration 
command  (e.g.  throttle  movement).  This  acceleration  command  is  assumed  to 
be  proportional  to  airspeed  errors  of  20  knots  or  less;  for  error  magnitudes  of  more 
than  20  knots,  maximum  acceleration  (ACR)  or  deceleration  (DCR)  for  the  aircraft 
type  is  commanded. 

Due  to  pilot-aircraft-engine  response  time,  the  commanded  acceleration/ 

deceleration  is  not  immediately  achieved.  The  inner  loop  from  A to  A represents 

c 
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a first-order  lag  response  with  a time  constant  of  t^,  . The  block  labeled  "Energy 
Constraint"  limits  the  actual  longitudinal  acceleration  (or  deceleration)  if  the 
aircraft  is  simultaneously  climbing  (or  descending).  The  acceleration  limit  is 
determined  by  the  excess  of  the  maximum  rate  of  climb  (or  descent)  capability 
over  the  actual  rate.  When  the  aircraft  is  climbing  at  its  maximum  capability, 
it  has  no  additional  power  available  for  acceleration.  Similarly,  during  a maximum 
rate  of  descent,  the  pilot  cannot  simultaneously  reduce  airspeed. 

The  next  block  in  Figure  6 resolves  the  longitudinal  and  lateral  accel- 
erations into  components  along  and  normal  to  the  ground  referenced  velocity. 

These  are  integrated  to  obtain  V and  X,  which  are  used  with  the  existing  wind 

9 

magnitude  and  direction  to  give  the  true  airspeed  and  heading  angle.  In  general, 
the  pilot  does  not  control  true  airspeed,  but  indicated  airspeed,  which  is  based 
on  the  difference  between  dynamic  and  static  pressure  measured  by  the  air  data 
system . 

The  groundspeed  is  resolved  into  North  and  East  components,  which  are 
integrated  to  obtain  the  aircraft  coordinates  x,y  in  the  ground  reference  frame. 

The  derivatives  x and  y,  as  well  as  all  other  derivatives  in  the  target  dynamics 
models  are  expressed  in  appropriate  units  per  second. 

The  velocity  and  heading  response  models  could  be  slightly  simplified 
by  neglecting  the  effect  of  the  aircraft  "crab  angle"  - X when  flying  through  a 
crosswind.  This  would  decouple  the  groundspeed  and  ground  track  angle  equations. 
However,  since  the  crab  angle  can  often  be  20°  or  more,  this  simplification  could 
result  in  sizable  errors,  particularly  when  turning  in  a wind  environment. 

The  integration  of  the  equations  in  the  velocity  and  heading  response 
loops  of  Figure  6 is  considerably  more  complex  than  in  the  case  of  the  altitude 
response.  Closed-form  solutions,  while  possible,  would  be  extremely  cumbersome. 
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A simple  numerical  integration  technique,  such  as  Runge-kutta,  should  provide 

sufficient  accuracy.  An  alternative  procedure  would  be  to  assume  that  the  signal 

is  constant  over  the  time  step  interval  at  a selected  point  in  each  loop  (e.g., 

A and  0 ).  Then  simpler  analytical  expressions  could  be  obtained  for  the  remainder 
Xc  C 

of  the  loop . 

3.3  MACH  NUMBER  COMMAND 

In  most  jet  aircraft,  the  pilot  will  often  be  controlling  Mach  number  instead 
of  airspeed,  particularly  during  the  climbout  and  cruise  phases.  The  control 
technique  is  essentially  the  same  as  for  airspeed  control,  with  the  exception 
that  the  longitudinal  acceleration  is  adjusted  according  to  the  error  in  Mach 
number.  Figure  7 shows  the  control  logic  for  the  Mach  number  command  loop. 


Figure  7.  Block  Diagram  of  Mach  Number  Command  Loop. 
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During  this  mode  of  operation,  the  A from  Figure  7 is  used  in  place  of  that 

c 

generated  in  Figure  6.  Otherwise,  the  control  loop  of  Figure  6 is  unchanged 
for  Mach  number  commands. 

The  Mach  number  tolerance  shown  is  0.005,  which  corresponds  to  a true 
airspeed  tolerance  of  the  order  of  3 knots.  Linear  control  is  assumed  for  Mach 
number  errors  between  0.005  and  0.10,  and  full  control  for  errors  in  excess  of 
0.10  Mach  number . 

3.4  HEADING  COMMAND 

The  heading  command  loop  shown  previously  in  Figure  6 applies  whenever 
the  pilot  is  attempting  to  maintain  a desired  magnetic  heading,  as  when  he  is 
being  vectored  by  ATC  for  example. 

Any  difference  between  desired  and  indicated  heading  is  used  to  generate 
a bank  angle  command  0 , according  to  the  block  labeled  "bank  angle  logic". 

A very  close  tolerance  on  heading  is  usually  maintained;  the  value  selected  in 
Figure  6 is  shown  as  0.5°.  Any  heading  errors  greater  than  this  are  corrected. 

The  maximum  commanded  bank  angle  a depends  upon  the  size  of  the  initial  heading 

error,  e,(0),  and  the  standard-rate  turn  bank  angle,  tan  \v  ai/g).  If  the  desired 

\|i  a 

heading  change  is  less  than  the  standard-rate  turn  bank  angle,  the  maximum  bank 
angle  command  is  equal  to  the  initial  heading  error;  otherwise,  the  bank  angle 
for  a standard-rate  turn  is  used.  The  initial  bank  angle  command  is  maintained 
until  the  heading  error  is  reduced  to  one-half  the  bank  angle  command  a/2  (Ref.  13). 
Then  linear  bank  control  is  used  until  the  heading  error  is  within  tolerance.  The 
value  of  a should  be  recomputed  whenever  a new  heading  command  is  issued. 

The  actual  bank  angle  is  approximated  by  a first-order  lag  response  to 
the  commanded  value.  The  lag  time  constant  approximates  the  pilot  and  aircraft 
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response  to  the  bank  command.  The  resulting  bank  angle  produces  a lateral  ac- 
celeration by  rotating  the  lift  vector  out  of  the  vertical  plane.  If  the  aircraft 
ang le-of-attack  is  adjusted  to  keep  the  vertical  component  of  lift  equal  to  the 
aircraft  weight,  the  lateral  acceleration  is  = g tan  0. 

As  explained  previously,  the  longitudinal  and  lateral  accelerations 

are  resolved  into  ground  referenced  coordinates  and  integrated  to  obtain  V and  X, 

9 

which  then  determine  V and  ty.  NOTE:  ground  track  and  heading  angles  must 

-t 

be  restricted  to  the  range 

0°  < X , \|[  < 360° 

However,  in  the  integration  of  the  X equation,  it  may  be  necessary  to  use  an 
intermediate  variable  which  is  unrestricted,  and  then  add  or  subtract  360°  as 
required  to  obtain  X. 

3.5  COURSE  TRACKING  COMMAND 

During  the  greatest  part  of  most  flights,  it  is  necessary  for  the  pilot  to 

follow  a given  course  over  the  ground  by  reference  to  various  navigation  systems. 

The  most  frequently-used  navaid  for  at  least  the  near  future  is  the  VOR/DME 

system  which  defines  the  low  altitude  Victor  Airways  and  high  altitude  Jet  Routes. 

Figure  8 defines  the  nomenclature  used  for  the  course  tracking  model. 

The  figure  shows  the  situation  of  an  aircraft  attempting  to  track  outbound  on  the 

9^  radial  of  a VOR/DME  station.  The  object  of  the  course  tracking  loop  is  to 

null  the  cross  track  distance  d = R sin  (9  - 9).  However,  to  prevent  large  over- 

c 

shoots  of  the  desired  track,  some  rate  damping  must  be  provided,  either  in  the  form 
of  d or  ground  track  error,  X - 0 . Since  in  most  cases  neither  of  these  is  readily 
available  aboard  the  aircraft,  the  difference  between  heading  and  desired  track. 
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Figure  8.  Geometry  and  Nomenclature  for  Course  Tracking  Model. 

\|i  -0  (i.e.,  the  zero-wind  ground  track  error),  is  frequently  used  as  a substitute. 

A model  for  this  process  is  analyzed  in  Kayton  and  Fried  (Ref.  14) 
and  is  shown  in  block  diagram  form  in  Figure  9.  The  analysis  of  this  model  neglects 
the  bank  angle  delay,  which  is  of  the  order  of  1 - 2 seconds,  and  assumes  the 
aircraft  is  in  the  vicinity  of  the  desired  course.  The  equations  of  motion  for  the 
no-wind  response  of  the  tracking  loop  are  linearized,  to  produce  a second-order 
system  with  natural  frequency  uu ^ and  damping  ratio  The  values  of  and 
can  be  adjusted  by  the  choice  of  the  gains  K-j  and  as  indicated  in  the  figure. 

A reasonable  response  is  obtained  with  uo^  —2  rad/min  and  - 0.7,  which  is 
equivalent  to  setting  the  gains 
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Figure  9.  Block  Diagram  of  Course  Tracking  Model. 
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NAVIGATION 

SYSTEM  STANDARD  -RATE 

TURN  LIMIT 


K j -12.0  deg/nm 

K0  = 0.000204  V nm/deg 

l a 

The  effect  of  the  wind  wil!  be  to  cause  slight  variations  typical  of  the 

variation  between  pilots  or  autopilots.  The  value  of  K_  is  made  a function  of  V 

2 a 

to  avoid  specifying  different  gains  for  various  speed  aircraft.  The  commanded 
bank  angle  for  the  course  tracking  model  must  be  limited  to  prevent  exceeding 
a standard  rate  turn,  where  a Is  the  magnitude  of  the  bank  angle  that  produces 
a standard  rate  turn.  The  resulting  bank  angle  command  0^  is  then  used  in  the 
heading  loop  of  Figure  6 in  place  of  the  value  shown  there  for  heading  commands. 

-i 

The  dynamics  of  the  bank  angle  response  introduce  additional  variations  between 
aircraft  which  are  typical  of  the  real  world.  As  a result  the  actual  response  will 
differ  slightly  from  that  of  a linear,  second-order  system  with  natural  frequency 
uu  and  damping  ratio 

The  model  for  tracking  with  a navigation  system  that  displays  off-track 
displacement  as  a linear  distance  instead  of  an  angle  is  only  slightly  different  from 
the  course  tracking  model  discussed  above.  This  could  be  used,  for  example,  with 
an  inertial  navigator  or  an  RNAV  system.  The  cross-track  error,  d,  would 
come  directly  from  the  navigation  system.  The  damping  signal  could  also  be 
obtained  directly  from  the  navigation  system  but  that  is  not  necessary. 

The  tracking  model  in  Figure  9 will  settle  down  to  zero  cross  track  error 
if  there  are  no  winds.  However,  in  the  presence  of  a cross  wind,  the  steady-state 
crab  angle  will  result  in  a small  cross  track  bias  error  d^  = - K 2(Xss  ' ♦„)•  This 
bias  could  be  eliminated  by  means  of  a washout  filter  or  by  adding  a term  proportional 
to  the  integral  of  d to  the  bank  angle  command.  However,  the  additional  complexity 
does  not  seem  to  be  warranted  for  the  next  stage  in  improving  the  ATC  simulation. 
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3.6 


COURSE  INTERCEPT  COMMAND 


In  most  instances,  the  aircraft  will  be  either  vectored  to  the  desired 
course,  or  transferring  from  one  leg  of  a route  to  another  with  relatively  small 
changes  in  track  angle.  In  these  situations,  the  tracking  model  will  be  perfectly 
satisfactory.  However,  the  pilot  will  occassional ly  be  directed  to  intercept  and 
maintain  a given  course  which  may  be  outside  the  valid  region  of  the  linear  tracking 
model.  In  the  latter  situations,  some  logic  is  required  to  guarantee  proper  intercept 
of  the  commanded  course. 

In  most  course  intercept  situations.  It  is  expedient  to  achieve  track  with 

the  least  delay.  For  the  idealized  case  of  constant  airspeed,  constant  altitude, 

no  winds  and  instantaneous  bank  angle  response,  the  minimum-time  intercept  paths 

consist  of  segments  of  circles  and  straight  lines  (Ref.  15).  A closed-loop  feedback 

control  chart  can  be  obtained  for  this  case  showing  the  optimum  bank  angle 

(+  o,0,-  a)  as  a function  of  cross  track  error,  d,  and  ground  track  error,  0 - \|i 

c 

(see  Fig.  8).  In  the  more  realistic  situation  being  considered  here,  a simpler 

sub-optimum  control  chart  can  also  be  defined  which  will  always  place  the  aircraft 

near  the  desired  course  where  the  course  tracking  model  can  be  correctly  applied. 

The  sub-optimum  minimum-time  course  intercept  logic  is  shown  in  Figure  10. 

The  ordinate  has  been  normalized  with  respect  to  the  radius  of  a standard-rate 

turn,  V /uu . In  the  regions  labeled  "track  mode"  the  course  tracking  model  discussed 

above  will  be  used  to  determine  the  bank  angle  command.  Note  these  regions 

extend  to  two  turn  radii  either  side  of  the  desired  course,  and  to  a heading  angle 

wi  thin  90°  of  the  course.  Outside  of  these  regions,  full  right  bank  (0^  - + a), 

full  left  bank  (0  = - a)  or  straight  flight  (0  = 0)  is  commanded  until  the  track 

c c 

mode  regions  are  penetrated.  For  large  initial  cross  track  errors,  this  logic  will 
call  for  a standard  rate  turn  of  minimum  heading  change  until  the  aircraft  is  headed 
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Figure  10.  Bank  Angle  Logic  for  Course  Infercept. 


nearly  perpendicular  to  the  course;  the  aircraft  will  then  fly  straight  until  it  is 
within  2 turn  radii,  after  which  it  will  switch  into  the  tracking  mode  of  operation. 

The  width  of  the  0c  = 0 regions  in  Figure  10  is  adjustable,  but  10°  seems 
to  be  satisfactory.  These  regions  are  expansions  of  the  straight  lines  at  0^  - \|i  = 90° 
and  270°  in  the  idealized  optimum  situation.  The  size  of  the  track  mode  region 
was  somewhat  arbitrarily  chosen  for  ease  of  implementation.  In  actuality,  maximum 
bank  angle  will  be  commanded  in  a large  portion  of  this  region  due  to  the  bank 
angle  constraints  in  Figure  9. 

Three  example  course  intercept  trajectories  are  shown  in  Figure  10  to 
illustrate  the  performance  of  the  logic.  The  ground  tracks  for  these  trajectories  ' 
relative  to  the  desired  course  are  sketched  in  the  upper  portion  of  the  figure. 

For  comparison,  the  corresponding  optimum  turns  are  also  sketched.  The  major 
differences  between  the  two  paths  are  a delay  in  achieving  the  standard  turn 
rate  for  the  suboptimum  case,  and  a lower  frequency  response  with  a slight  overshoot 
within  the  tracking  mode  region. 

3.7  NOMINAL  TARGET  PARAMETERS 

A number  of  aircraft-related  parameters  are  required  to  utilize  the  target 
dynamics  model  just  described.  Appropriate  values  of  these  various  parameters 
should  be  specified  by  the  simulation  user  for  his  particular  investigation.  However 
as  a guide  to  the  user,  Table  1 contains  suggested  nominal  values  of  the  parameters 
for  a variety  of  aircraft  types.  As  explained  earlier  in  this  section,  the  parameters 
in  Table  1 illustrate  typical  estimated  values  for  normal  enroute  operations.  These 
values  should  not  be  utilized  in  general  for  relatively  quick  reaction  time  situations 
such  as  collision  avoidance.  To  the  extent  possible,  the  aircraft  types  and  parameters 
presented  are  consistent  with  those  used  in  previous  DOT  simulations  at  TSC  and 
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Table  1 . Typical  Aircraft  Parameters  for  Target  Dynamics  Model 
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NAFEC.  These  values  could  be  permanently  included  in  the  simulation  to  provide 
default  values  for  those  quantities  not  specified  by  the  user. 

As  mentioned  earlier  no  two  aircraft,  even  of  the  same  type,  are  truly 
identical  because  of  condition,  equipment,  loading  and  other  factors.  Consequently, 
the  parameters  presented  in  Table  1 should  be  considered  as  average  values.  The 
specific  values  for  a particular  simulation  target  should  be  selected  from  a normal 
random  distribution,  having  a mean  as  shown  in  Table  1,  with  a suggested  standard 
deviation  of  10%  to  20%  of  the  mean  value.  However,  care  must  be  taken  that 
the  selected  parameters  are  always  positive. 
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SECTION  4 


AIR  DATA  SYSTEM  MODEL 

The  aircraft  air  data  system  consists  of  aerodynamic  and  thermodynamic 
sensors,  some  form  of  computer  and  various  displays.  The  outputs  of  the  air  data 
system  are  indicated  values  of  various  flight  parameters  such  as  altitude,  air- 
speed or  Mach  number.  Although  it  is  not  generally  considered  a part  of  the  air 
data  system,  the  heading  indicator  will  be  included  in  the  air  data  system  model 
for  convenience . 

A general  block  diagram  of  the  air  data  system  is  shown  in  Figure  1 1 . 

A detailed  explanation  of  the  air  data  system  operation  is  not  necessary  for  the 
purposes  of  this  model;  such  a description  may  be  found  in  other  sources  (e.g.  Ref.  14) 
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Figure  1 1 . General  Block  Diagram  of  Air  Data  System  Model . 
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The  air  data  system  errors  generally  consist  of  scale  factor  errors  and  biases. 

Each  element  of  the  air  data  system  model  includes  a reasonably  simple  represen- 
tation of  the  overall  error  in  the  indicated  quantity. 

4.1  ALTIMETER 

The  barometric  altimeter  indicates  the  pressure  altitude  of  the  aircraft  based 
on  the  static  pressure  measurement  and  a standard  temperature-altitude  profile  of  the 
atmosphere.  Local  barometric  fluctuations  can  be  corrected  for  by  adjusting  the  sea- 
level  barometric  pressure  in  flight.  The  principal  altimeter  errors  are  due  to  static 
defect,  nonstandard  temperature-lapse  rate  and  instrument  errors. 

The  static  defect  is  the  difference  between  the  free-stream  static  pressure  and 
the  measured  value,  and  is  usually  not  significant  at  speeds  below  about  Mach  0.4. 

A static  defect  is  usually  calibrated  as  a function  of  Mach  number,  and  automatically 
corrected  for  by  the  air  data  system  (c.f.  Ref.  16).  A typical  maximum  static 
defect  is  350  ft,  with  a spread  of  100  ft  among  aircraft  of  the  same  type  (Ref.  14). 

A non-standard  temperature-lapse  will  result  in  an  error  between  true  altitude 
and  pressure  altitude,  but  all  aircraft  will  be  similarly  affected.  Therefore  lapse- 
rate  errors  will  not  appreciably  affect  the  relative  vertical  spacing  of  aircraft, 
although  terrain  clearance  can  be  in  error  as  much  as  several  tens  of  feet  per  thousand 
feet  of  true  altitude. 

Altimeter  instrument  errors,  including  transducer,  computer  and  indicator 
errors,  have  been  widely  discussed  (Refs.  17  - 20),  especially  for  altitudes  above 
30,000  ft.  The  altitude  variation  of  the  instrument  error  is  inversely  proportional 
to  the  atmospheric  density.  The  SAE  design  objective  is  an  accuracy  of  0.2 
percent  or  25  ft,  whichever  is  greater  (Ref.  21).  The  latest  airline  specifications 
(ARINC  575-3)  call  for  a 3 a instrument  error  ranging  from  15  ft  at  sea  level  to 
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80  ft  at  50,000  ft  (Ref.  22).  It  is  estimated  that  the  standard  deviation  of  all 
present  altimeters  is  about  25  ft  at  sea  level . 

The  altimeter  model  is  based  principally  on  the  instrument  error  discussed 
above.  The  indicated  altitude  is 

h|  = b + eh(h)  (9) 

where  the  error  e^(h)  is  given  by 

*h(h)=  €h(0)[p(h)/p0]_1  OO) 

The  atmospheric  density  ratio  can  be  closely  approximated  as  an  experimental  function 
of  altitude: 

-h/h, 

p(h)/  Pq  = e ^ 1 ^ 

where  h^  =31,186  ft.  The  sea  level  instrument  error  e^(Q)  for  a particular  altimeter 
should  be  selected  from  a norma!  distribution  with  zero  mean  and  a standard  deviation 
of  about  25  ft: 

eh(0)  = N[0,25]  0?) 

4.2  VERTICAL  SPEED  INDICATOR 

The  standard  vertical  speed  indicator  (VSI),  or  rate  of  climb  indicator,  senses 

differential  changes  in  the  static  pressure  by  means  of  a calibrated  leak.  As  static 
pressure  changes  during  entry  to  a climb  or  descent,  the  VSI  immediately  shows  a 

change  of  vertical  direction.  However,  until  the  differential  pressure  stabilizes, 
reliable  rate  indications  cannot  be  obtained.  A 6-  to  9-second  lag  normally  is 
required  to  equalize  or  stabilize  the  pressures.  The  instantaneous  vertical  speed 
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indicator  is  a recent  development  which  incorporates  acceleration  pumps  to  stabilize 
the  pressure  differential  without  the  usual  lag  time.  Typical  VSI  errors  are  estimated 
to  be  about  5%  to  10%  of  the  altitude  rate.  ARINC  Specification  575-3  (Ref.  22) 
calls  for  an  accuracy  of  30  ft/min  or  5%,  whichever  is  greater. 

The  uncertainty  in  a particular  VSI  is  modeled  as  a fixed  scale  factor 

error: 

h,  = h+€^=K^h  (13) 

where  the  scale  factor  is  selected  from  a normal  distribution  with  a a of  5%: 

= N[1  .0,0.05]  (14) 

Thus,  the  absolute  error  in  the  vertical  speed  indicator  is  given  by 

eh  = (Kh  “ 1)^  (15) 

The  pressure  stabilization  lag  time  is  not  included  in  the  model,  since  altitude 
rate  is  not  controlled  very  tightly  during  transient  conditions,  and  the  instantaneous 
VSI  is  available  if  required. 


4.3  AIRSPEED  INDICATOR 

True  airspeed,  V , is  the  actual  speed  of  the  aircraft  relative  to  the 
surrounding  air  mass.  Of  more  importance  to  the  pilot,  however,  is  indicated 
airspeed,  V|^<-,  which  determines  maximum  and  minimum  safe  flying  speeds, 
and  on  which  climb,  descent  and  cruising  speed  schedules  are  often  based.  V|^ 
is  measured  directly  by  the  difference  between  the  impact  pressure,  p^.,  sensed 
by  the  pitot  tube,  and  the  static  pressure,  ps: 


V |AS  = 0.5925  ^2(pt  - ps)/p0  (kts) 


(16) 
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where  pQ  is  the  standard  density  at  sea  level.  Indicated  airspeed  is  related  to 
true  airspeed  as  follows.  The  static  defect  and  instrument  errors  correct  to 

calibrated  airspeed  $ is  multiplied  by  ^Pq/p  to  correct  for  the  density 

at  the  aircraft  altitude,  resulting  in  equivalent  airspeed.  A correction  which 
accounts  for  the  compressibility  of  the  air  at  high  altitude-speed  regimes  corrects 
equivalent  airspeed  to  true  airspeed,  V . Below  Mach  0.8  at  30,000  ft,  equivalent 
airspeed  and  true  airspeed  are  identical  within  7 percent  (Ref.  14). 

The  model  for  obtaining  indicated  airspeed  from  true  airspeed  is  shown 
in  Figure  12.  The  compressibility  correction  is  neglected  since  It  is  small  relative 


Figure  12.  Block  Diagram  of  Indicated  Airspeed  Model. 


to  the  density  correction,  and  since  Mach  number  rather  than  V|^<^  is  normally 
controlled  whenever  the  compressibility  is  significant.  Consequently,  true  air- 
speed is  converted  to  by  fbe  density  ratio,  which  can  be  closely  approximated 

as  shown  in  Equation  (11).  The  instrument  error  which  distinguishes  between 
and  is  approximately  proportional  to  airspeed.  Thus,  for  a particular 

aircraft,  the  error  in  can  be  represented  by  a scale  factor  Ky  selected 

from  a normal  distribution: 
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Ky  = N[  1 .0,0.02] 


07) 


where  the  standard  deviation  of  2%  is  estimated  to  be  typical. 

4.4  MACH  INDICATOR 

Mach  number  is  required  in  high  performance  aircraft  for  defining  safe 
speed  ranges,  fuel  management,  jet-engine  control,  and  calculating  other  aero- 
dynamic quantities.  Mach  number  is  the  ratio  of  the  aircraft's  true  airspeed  to 
the  local  speed  of  sound.  It  is  generally  calculated  in  flight  from  the  ratio  of 
impact  to  static  air  pressure  (Ref.  14). 

Figure  13  shows  the  model  for  indicated  Mach  number,  M| . The  local 


M 


Figure  13.  Block  Diagram  of  Indicated  Mach  Number  Model. 

speed  of  sound,  V^h)  in  knots,  can  be  closely  approximated  by  the  expression 

( 1 .734472  a/i45,501  - h , h<  36,089  ft 
V (h)  = (18) 

s (573.8  , h>  36,089  ft 

The  instrument  error  in  the  Mach  indicator  is  approximately  2%,  which  can  be 
represented  as  another  scale  factor. 
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Kk.  = N[  1 .0,0.02] 

M 


09) 


4.5  HEADING  INDICATOR 

Although  the  heading  indicator  is  not  truly  a part  of  the  air  data  system, 
it  is  included  here  for  convenience.  The  gyroscopic  heading  indicator  provides 
the  pilot's  primary  direction  information.  In  most  high  performance  aircraft  the 
heading  indicator  is  slaved  to  a remote  compass  which  continuously  maintains 
the  north  reference.  Less  expensive  gyros,  used  in  most  general  aviation  aircraft, 
must  be  manually  reset  from  the  magnetic  compass  periodically.  Thus,  the  average 
accuracy  of  the  heading  indicators  varies  considerably  among  various  aircraft. 

The  principal  errors  are  due  to  inaccuracies  in  the  magnetic  compass  reference 
and,  for  non-slaved  instruments,  to  gyro  drift  between  resettings. 

The  indicated  heading  model  assumes  a constant  error  selected  from 


a normal  distribution  with  zero  mean: 


♦ , = ♦ + 6 


* 


% = n[°,g^1 


(20) 

(21) 


The  standard  deviation  of  the  heading  error,  a , is  estimated  to  range  between 
about  1°  (for  commercial  jets)  to  3°  or  more  (for  smaller  general  aviation  aircraft). 
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SECTION  5 


NAVIGATION  SYSTEMS  MODEL 

A great  number  of  navigation  systems  have  been  developed  to  provide 
the  pilot  with  position  information  and/or  course  guidance  (cf.  Ref.  14).  These 
range  from  ground-based  radio  direction  finding  to  sophisticated  self-contained 
inertial  navigation  systems.  By  far  the  most  widely  used  system  at  present,  and 
for  at  least  the  near  future,  is  the  VOR/DME  system,  upon  which  the  US  low- 
altitude  Victor  Airways  and  high-altitude  Jet  Routes  are  based.  An  increasing 
number  of  airborne  area  navigation  (RNAV)  systems  are  coming  into  operation  to 
provide  much  greater  flexibility  from  the  existing  VOR/DME  system. 

5.1  VOR/DME  SYSTEM 

This  system  is  comprixed  of  two  separate  but  complementary  radio  navaids 
which  are  normally  used  in  conjunction  with  one  another: 

• Very-high-frequency  Omnidirectional  Range  (VOR) 

• Distance  Measuring  Equipment  (DME) 

The  development  and  operation  of  these  systems  has  been  widely  documented  in 
the  literature  (e.g.  Refs.  23  - 30)  and  will  not  be  repeated  here.  The  VOR 
indicates  the  magnetic  bearing  of  the  aircraft  from  the  selected  ground  station, 
while  the  DME  measures  the  slant  range  between  them  (see  Fig.  14). 

VOR  Model 

The  correct  magnetic  bearing  from  the  VOR  station  to  the  aircraft  is: 

e = + tan_1  " \)/(y„  “ ys)]  ^ 

M a s a s 


- 39  - 


TRUE 


Figure  14.  Geometry  and  Nomenclature  forVOR/DME  Model. 


The  indicated  VOR  bearing  is  in  error  by  the  amount  e^:  ' 

V 


6y  — 9 + €x 


(23) 


where  the  total  error  includes  receiver  and  transmitter  biases  and  a course  roughness 

error: 


ev-  cR+  eT+  ec 


(24) 


For  each  aircraft,  is  a random  constant,  selected  from  a normal 
distribution  with  zero  mean  and  a standard  deviation  a^,  i.e.: 

= N[0,a^] 

Values  of  for  high  quality  receivers  are  approximately  0.3°  (Ref.  28)  to  0.6° 


(25) 
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(Ref.  30)  with  0.2°  claimed  as  possible  for  future  ARINC  equipment  (Ref.  31). 

Typical  values  of  for  average  receivers  seem  to  be  about  1.0°  to  1 .5°  (Refs. 

10,  14,  24,  27). 

Is  also  a random  constant  for  each  VOR  station: 

ey  = N[  0,0-j.  ] (26) 

For  conventional  VOR  (CVOR)  stations,  Oy  is  typically  around  0.35°  to  1.0° 

(Refs.  10,  14,  23,  24,  28  - 30).  For  doppler  VOR  (DVOR)  transmitters  is 
about  0.05°  to  0.25°  (Refs.  14,  24,  27,  32).  Suggested  values  are  0.9°  for 
CVOR  and  0.15°  for  DVOR. 

The  position-dependent  course  roughness  error  is  modeled  as  colored  noise 
whose  correlation  time  is  inversely  proportional  to  aircraft  groundspeed,  i.e. 

9 -O)  V | T | 

<£*[e  (t)e  (t  + t)]  = a e c 9 (27) 

CO  c 

The  standard  deviation  a ^ is  approximately  0.9°  for  CVOR,  0.35°  for  DVOR 

and  0.2°  for  precision  DVOR  (Refs.  24,  28,  32,  33).  The  course  roughness 

-3 

frequency  is  estimated  between  0.2  and  2.0  X 10  rad/sec/kt  for  CVOR  and 
about  4x10^  rad/sec/kt  for  DVOR  (Ref.  28).  A suggested  value  for  CVOR  is 
1 X 10  2 rad/sec/kt. 

The  course  roughness  error  is  generated  from  Gaussian  white  noise  as 
shown  in  Figure  15.  The  autocorrelation  function  of  the  driving  white  noise  is 

<?[wy(t)wy(t  + t)]  = (2ac2/ajcvg)f>(T)  (28) 

The  VOR  model  could  be  simplified  by  ignoring  the  magnetic  variation 

8 and  usinq  true  bearing  as  the  measured  quantity.  If  this  is  done,  true  heading 
M 

instead  of  magnetic  heading  must  be  used  for  tracking  a radial  in  the  target  dynamics 
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Figure  15.  Block  Diagram  of  VOR  Error  Model. 


model.  Thus  it  would  not  be  necessary  to  maintain  a continuous  indication  of 
6^  at  each  aircraft's  position. 

The  frequency  of  the  VOR  course  roughness  error  model  is  independent 
of  aircraft  position  or  heading.  This  is  not  strictly  accurate,  but  it  does  provide 
a substantially  simpler  model.  In  actuality,  the  course  roughness  fora  specific 
VOR  is  definitely  position-related,  since  much  of  it  is  the  result  of  signal  attenuation 
and/or  reflection  by  the  local  terrain,  structures  and  nearby  aircraft.  For  example, 
two  aircraft  following  nearly  the  same  path  within  a short  period  of  time  would 
probably  experience  very  similar  course  roughness.  Also,  the  course  roughness 
frequency  would  generally  be  different  when  flying  along  radials  and  while  orbiting 
the  station,  as  well  as  vary  with  altitude.  However,  due  to  the  limited  amount 
of  data  and  the  uncertainty  associated  with  it,  the  simpler  model  was  chosen. 

A much  more  complex  model  and  a large  amount  of  empirical  data  would  be  necessary 
to  simulate  the  range,  bearing  and  altitude  dependence  of  the  course  roughness 
at  a particular  VOR  station. 
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DME  Model 


The  true  slant  range  in  nm  between  the  DME  station  and  the  aircraft  is 


R = ^xa  “ xs)2  + (yQ  - ys)2  + [0.0001646(hQ  - hs)]2  (29) 

where  x and  y are  in  nm  and  h is  in  ft.  The  slant  range  indicated  by  the  DME 
equipment  is  in  error  by 

Rn  = R + en  (30) 


The  slowly-varying  DME  range  error  is  modelled  as  colored  noise  with  a correlation 
time  of  several  minutes: 


^eD^eD^  + = 


O. 


2 -HA 


D 


(31) 


Tp  estimates  range  between  300  sec  and  600  sec  (Refs.  25,34);  400  sec  is  recommended 
as  a compromise  value. 

The  standard  deviation  in  includes  independent  contributions  from  the  air- 
craft interrogator  (oj)  and  the  DME  station  transponder  ( a <-): 


°D 


(32) 


a | estimates  range  from  less  than  0.1  nm  to  about  0.25  nm,  while  a<-  is  typically 
between  0.03  nm  and  0.1  nm  (Refs.  14,  24,  28,  30,  35,  36).  The  overall  DME 
error  o ^ estimates  range  from  0.015  nm  to  0.5  nm  (Refs.  10,  26,  34,  37).  The 
airlines'  goal  for  is  ± 83  ft  (.014  nm)  for  all-weather  landing  operations  (Ref.  38). 
The  recommended  value  of  for  most  operators  is  0.1  nm. 

As  shown  in  Figure  16,  the  DME  correlated  error  is  generated  from  white 
noise  with  the  autocorrelation  function 


- 43  - 


R (t) 


£[wD(t)wD(t  + t)]  = 2aDTd6(T)  (33) 

Because  the  DME  error  is  colored  noise  with  a long  correlation  time,  it  could 
be  closely  approximated  as  a constant  bias  rather  than  a slowly-varying  bias.  This 
would  substantially  reduce  the  computation  time  for  the  model.  Further  justifications 
for  this  simplification  are  that  the  DME  error  is  generally  much  smaller  than  the  VOR 
error,  and  that  very  little  information  exists  on  the  correlation  time  of  the  error. 


Limits  (Fig . 17) 

The  VOR  and  DME  signals  are  both  limited  to  I ine-of-sight  to  the  station. 
Assuming  a circular  earth,  this  restricts  the  aircraft  altitude  to  approximately: 


h > 

a 


0.788  ^(xa  - »s)2  + (ya  - ys)2  - Vh 


(34) 


For  adequate  VOR  or  DME  signals,  the  slant  range  must  be  within  a maximum 


value,  R<  R . For  VOR's,  R is  determined  to  insure  freedom  of  interference  ‘ 
max  max 
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AIRCRAFT 


by  another  station  transmitting  on  the  same  frequency.  For  high  altitude  class 
VORs,  R x is  130  nm  (Ref.  39).  The  maximum  ranges  of  most  DMEs  vary  from 
100  to  300  nm  (Refs.  14,  26).  ARSNC  568  (Ref.  38)  specifies  a maximum  error 
of  0.2  nm  at  a range  of  380  nm  for  future  airline  equipment. 

For  a usable  VOR  signal,  the  aircraft  must  be  outside  the  cone  of  ambiguity 
over  the  station.  This  limits  the  altitude  to 

h < h + 6,076. 1 \(x  - x )^  + (y  - y )^  cot  y (35) 

as  las  'a  's 

where  the  half-cone  angle  y is  typically  between  15  and  30  (Refs.  23,  28,  31). 

An  alternative  to  explicitly  including  the  limits  in  the  model  is  for  the 
model  user  to  include  them  implicitly.  Under  this  approach,  the  user  could  apply 
the  model  within  the  region  of  validity,  or  would  use  the  limitations  to  interpret 
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his  results.  This  could  reduce  the  computation  time  slightly  by  avoiding  the  three 
checks  every  update. 

5.2  AREA  NAVIGATION  SYSTEMS 

Although  the  VOR/DME  system  provides  nearly  complete  coverage  of 
US  airspace,  available  routes  are  limited  by  the  requirement  to  fly  from  station 
to  station.  This  has  several  disadvantages:  the  indirect  routes  are  longer,  much 
of  the  airspace  is  wasted,  and  the  danger  of  collision  is  increased  near  the  station 
where  numerous  airways  converge.  The  use  of  area  navigation  (RNAV)  equipment 
can  alleviate  these  problems  by  permitting  direct  point-to-point  navigation. 

RNAV  devices  have  been  under  development  for  many  years,  but  only 
within  the  past  half-decade  has  serious  testing  demonstrated  the  feasibility  of 
area  navigation  (Refs.  25,  40).  Since  then,  considerable  progress  has  been  made 
not  only  in  the  equipment  itself,  but  in  the  development  of  standards  and  procedures 
for  the  implementation  of  RNAV  in  the  National  Airspace  System  (Refs.  41  - 45) 

The  RNAV  system  computer  uses  standard  VOR/DME  information  to  determine 
the  aircraft's  position  relative  to  any  course  selected  by  the  pilot.  There  are  two  forms 
of  output.  The  more  common  output  is  range  and  bearing  to  a waypoint  which  has  been 
established  by  a radial  and  distance  from  a given  VOR/DME  station.  The  second  form 
of  output  is  a linear  display  which  shows  the  aircraft  position  relative  to  some  map 
coordinates.  The  former  approach  has  been  selected  for  the  VOR/DME  RNAV  model, 
while  the  latter  will  be  used  for  the  DME/DME  RNAV  model. 

VOR/DME  RNAV  Model 

Most  available  RNAV  systems  use  the  information  from  a single  VOR/DME 
station  for  navigation.  The  RNAV  routes  are  normally  specified  by  waypoints,  which 
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are  defined  by  their  bearing  and  distance  from  a given  VOR/DME  station.  These 
way-points  are  treated  as  conveniently-located  phantom  VOR/DME  stations,  and 
then  used  for  navigation  in  the  normal  fashion. 

Figure  18  illustrates  the  nomenclature  and  geometry. used  for  the  VOR/DME 


Figure  18.  Geometry  and  Nomenclature  for  VOR/DME  RNAV  Model. 


RNAV  model.  The  DME  indicated  slant  range  from  the  station  to  the  aircraft  is 


corrected  for  altitude 
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(36) 


d = y Rq2  * 10.0001646(ha  - hs)]2 

Slant  range  correction  is  norma!  but  not  universal  in  existing  RNAV  equipment.  The 
FAA  intends  future  rule-making  to  require  it.  Until  then,  RNAV  airways  will  maintain 
sufficient  separation  from  VOR/DME  stations  so  that  the  slant  range  correction 
is  negligible.  RNAV  approach  procedures  which  require  slant  range  correction 
will  be  limited  to  those  aircraft  which  are  so  equipped. 

The  distance  and  bearing  from  the  waypoint  to  the  aircraft  are  found 
geometrical  ly 


d ' = + d ^ _ 2dd  cos  (9  - 8,,) 

J W W W V 


9=9  - sin 

w 


. -1 


-Sin(0w  - V 


(37) 

(38) 


No  error  is  included  for  the  RNAV  computations  other  than  what  exists 
for  the  VOR/DME  sensors.  In  the  less  expensive  systems  there  might  be  some 
computational  error  to  be  added.  However,  the  more  sophisticated  RNAV  systems 
would  have  reduced  total  error  because  of  optimum  filtering  with  air  data  or  inertial 

measurements. 

If  the  VOR/DME  limits  are  included  in  the  model,  the  RNAV  inputs  may 
be  lost,  for  example,  when  flying  through  the  VOR  cone  of  ambiguity.  If  other 
navigation  information  is  not  provided  to  the  RNAV  system,  this  problem  can  be 
alleviated  by  continuously  estimating  groundspeed  from  position  differences  until 
the  signal  is  lost,  and  then  extrapolating  the  position  until  reacquisition  of  signal. 
Alternatively,  the  steering  logic  in  the  target  dynamics  model  would  maintain  a 
constant  heading  during  signal  loss. 


- 48  - 


DME/DME  RNAV  Model 


Instead  of  using  range  and  bearing  from  a single  VOR/DME  station,  the  RNAV 
computer  could  use  range  or  bearing  information  from  two  separate  stations.  Because 
the  DME  position  error  is  generally  much  smaller  than  the  VOR  position  error,  the  use 
of  two  DME's  could  theoretically  provide  much  higher  accuracy  than  two  VOR's  or  a 
single  VOR/DME. 

The  DME/DME  RNAV  model  described  below  does  not  represent  any  known 
existing  or  proposed  system.  However,  the  model  does  provide  an  indication  of  the 
possible  accuracy  of  a range-range-altitude  navigation  system.  The  DME/DME  RNAV 
model  geometry  is  illustrated  in  Figure  19.  The  two  indicated  DME  slant  ranges 
are  corrected  for  altitude 


INDICATED  AIRCRAFT  POSITION 


Figure  19.  Geometry  and  Nomenclature  for  DME/DME  RNAV  Model. 
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d = ^D1  - [0.0001646(hQ  -h^]2 


d2  = /RD2  " [°-0001646(ha  “ h2)]2 


The  baseline  horizontal  distance  between  the  two  DME  stations  is 


(39a) 


(39b) 


D= 


x^2  + (y2  - y^2 


(40) 


The  location  of  the  aircraft  is  determined  from  the  intersection  of  the  two  horizontal 
circles  of  radii  d-j  and  d 2*  centered  at  stations  1 and  2,  respectively.  Except  for  the 
special  case  d-j  + 62  < D,  there  are  two  possible  solutions  for  the  aircraft  position. 

To  resolve  this  ambiguity,  the  bearing  information  from  one  of  the  stations  (DME^) 
is  used . 

The  true  bearing  of  the  indicated  aircraft  position  from  the  VOR/DME^  station 
is  found  by  geometry: 


a = Oq  ± cos  ^ [(D2  + d^2  - d22)/2d^D] 
where  Oq  is  the  true  bearing  from  the  VOR/DME^  station  to  the  DME2  station: 


(41) 


<y, 


-1 


= tan"  (x2  " x1)/(y2  - y^ 


(42) 


Since  9y  - 0^  ~ a , the  choice  of  the  sign  in  Equation  (41)  is  determined  as  follows 

< Oq  — ♦ negative  sign 
<Oq  — . positive  sign 


9V  " 0M 


(43) 


Finally,  the  indicated  position  components  are  found: 


x = x,  + d,  sin  a 
a 1 1 


y = y , + d , cos  a 
' a '1  I 


(44a) 

(44b) 
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The  VOR  bearing  data  is  used  in  this  model  only  to  resolve  the  ambiguity 
of  position.  However,  the  redundant  information  provided  could  be  used  with 
an  optimal  filter  to  provide  somewhat  better  accuracy.  Likewise,  redundant  data 
from  other  navaids,  inertial  systems  and  the  air  data  system  can  be  filtered  together 
for  improved  performance  (Refs.  34,  35,  46,  47). 

Another  drawback  of  DME/DME  navigation  is  the  high  expense  of  the 
airborne  equipment,  in  comparison  to  VOR  equipment.  A single  airborne  set  could 
be  switched  between  two  or  more  ground  stations,  but  this  would  introduce  time 
delays  for  search  and  iock-on  of  each  station,  as  well  as  the  switching  complexity. 
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SECTION  6 


SURVEILLANCE  SYSTEM  MODEL 


The  present  ATC  system  uses  a dual  network  of  radar  systems  to  identify 
and  track  aircraft:  The  primary  system  uses  pulsed-type  search  radar;  the  secondary 
system  uses  pulsed  interrogators  and  airborne  transponders  to  enhance  the  return 
pulses  and  identify  the  targets.  Long  range  (200  nm)  air  route  surveillance  radars 
(ARSR)  supply  radar  data  for  en  route  control,  and  separate  shorter  range  (60  nm) 
airporK  survei  I lance  radars  (ASR)  supply  radar  data  in  the  critical  terminal  area. 
Table  2 summarizes  the  important  parameters  of  the  two  systems  (Ref.  10).  Computer 


Table  2.  Present  FAA  Radar  System  Parameters. 


Parameters 

Terminal  Radar 

En  Route  Radar 

Max  range 

60  nm 

200  nm 

Elevation  coverage 

0°  to  30° 

o 

o 

co 

o 

■4— 

o 

o 

Fan  beam  azimuth 
bandwidth 

1 .5° 

1 .35° 

Frequency 

2700  to  2900 
MHz 

1280  to  1350 
MHz 

Scan  Rate  (con- 
ti  nuous) 

15  rpm 

6 rpm 

PRF*  (nominal) 

1 200  pps 

360  pps 

Pulse  width 

0.833  |isec 

1 .95  psec 

i 

“ft 

Pulse  Repetition  Frequency 


processing  of  the  radar  surveillance  information  provides  easy  transfer  and  accurate 


processing  and  updating  of  flight  information,  aids  for  establishing  and  maintaining 
radar  identification  of  aircraft,  and  automatic  display  of  altitude  and  groundspeed 
information  with  aircraft  track  position. 

Control  of  IFR  traffic  is  based  primarily  on  use  of  the  ATC  radar  beacon 
system  (ATCRBS)  replies  for  surveillance.  Primary  radar  data  is  available  for 
tracking  unequipped  aircraft  and  for  airborne  transponder  failures.  Satellite  sur- 
veillance systems  are  presently  being  investigated  for  future  ATC  applications, 
particularly  over  oceanic  routes. 

6.1  GROUND  SURVEILLANCE 

The  ground  surveillance  model  described  below  is  based  on  the  ATCRBS 
in  the  Third  Generation  air  traffic  control  system.  The  geometry  and  nomenclature 
for  the  ground  surveillance  model  is  shown  In  Figure  20.  During  each  scan  the 


Figure  20.  Ground  Surveillance  Geometry  and  Nomenclature. 
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transponder  will  be  interrogated  several  times  as  the  radar  beam  sweeps  across 
the  aircraft.  The  beacon  processor  records  the  azimuth  and  range  of  each  successful 
beam  reply  (hit).  The  range  and  azimuth  of  each  hit  are  quantized  and  entered 
in  a two-dimensional  table  of  range-azimuth  cells,  as  illustrated  in  Figure  21. 
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Figure  21  . Range-Azimuth  Cells  in  Beacon  Target  Processor. 

A "1"  in  the  figure  indicates  a hit  was  detected  with  a range  between  R and  R + AR, 
and  with  an  azimuth  between  0 and  0 + A0.  In  the  terminal  area,  the  range  bin  aR 
is  1/16  nm,  and  the  azimuth  bin  A0  is  0.088°.  The  corresponding  values  for  the  en 
route  system  are  1/4  nm  and  0.225  . Based  on  certain  statistical  criteria,  a target 
detection  algorithm  scans  all  hits  at  a given  range  to  declare  a beacon  target  and 
to  establish  its  leading  and  trailing  edges  (Ref.  48).  The  azimuth  center  mark  of  the 
target  is  determined  by  averaging  the  leading  and  trailing  edge  positions  and  adding 
a bias  which  is  a function  of  the  target  detection  parameters. 

The  true  azimuth  from  the  radar  to  the  aircraft  is  (Fig.  19): 

0(t.)  = tan_1[(xa  - xr)/(yQ  - yf)]  (45) 

The  observation  times  t.  are  t.  = 0,T,2T,3T,  etc.  The  radar  scan  period  in 
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seconds  is 


T = 60/Q  (46) 

where  0.  is  the  radar  scan  rate  in  rpm. 

The  radar  azimuth  error  includes  the  effects  of  interrogator  and  trans- 
ponder errors  and  the  quantization  error  in  the  beacon  target  processor.  The  error 
is  modeled  by  assuming  the  unquantized  error  on  each  scan  is  a normally- 
distributed  random  variable  with  zero  mean  and  standard  deviation  a^: 

eX  = N[0'°A]  (47) 

a*  is  typically  0.25°  to  0.4°  for  the  terminal  radar,  with  0.12°  possible;  the 
slower  scan  rate  of  the  enroute  system  reduces  a ^ by  a factor  of  about  2.5  (Refs. 

10,  49,  50).  The  indicated  radar  azimuth  neglecting  the  quantization  error  is 
given  by 

e;(t5)  = e(t.)+  e^(t.)  (48) 

This  value  is  quantized  to  the  azimuth  bin  size  to  obtain  the  final  indicated  azimuth: 

9 (t.)  = n A9  (49a) 

where  n is  found  such  that 

n A0  < 9^(tj)  < (n  + 1 ) A9  (49b) 

The  true  range  in  nautical  miles  from  the  radar  to  the  aircraft  is 

R(t.)=W(x  -x)2  + (y  - y )2  + [0.0001646  (h  - h )2]  (50) 

where  x and  y are  in  nm  and  h is  in  ft.  The  indicated  radar  range  each  scan  is 
in  error  by  an  amount  due  to  the  quantization  effect  and  to  inaccuracies  in 
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the  elapsed  time  between  interrogation  and  reply.  The  primary  factor  in  the  range 

measurement  error  is  a bias  due  to  lack  of  precise  delay  control  in  the  aircraft 

transponder  system.  Therefore  the  range  measurement  error  neglecting  quantization 

is  modeled  as  a random  bias,  constant  for  each  aircraft,  which  is  selected 

from  a normal  distribution  with  zero  mean  and  standard  deviation  aD: 

R 

eR  = N[°,CTr]  (5,) 

Estimates  of  are  typically  1/32  to  1/16  nm  (190  ft  to  380  ft),  with  100  ft  possible 
with  tightened  equipment  specifications  (Refs.  10,  49,  50).  The  indicated  radar 
range  neglecting  the  quantization  error  is  given  by 

R;(t.)  = R(t.)  + eR  (52) 

This  value  is  quantized  to  the  range  bin  size,  to  obtain  the  final  indicated  value 

Rr(t.)  = m AR  (53a) 

where  m is  found  such  that 

m AR  < R^(t.)  < (m  + 1 ) AR  (53b) 

A possibly  significant  phenomenon  of  the  ATCRBS  is  the  probability  of 
lost  replies.  These  may  be  caused  by  1)  overinterrogation  by  many  interrogators 
which'cause  dead  time  in  the  transponders,  or  decrease  their  sensitivity;  2)  poor 
coverage  produced  by  interrogator  siting  limitations  and/or  multipath  nulls;  and 
3)  aircraft  antenna  shielding,  particularly  during  turning  maneuvers.  Dropout 
due  to  aircraft  banking  is  the  most  important  effect  for  the  simulation  since  the 
probability  of  a dropout  is  two  or  three  orders  of  magnitude  higher,  and  the  dropout 
is  likely  to  continue  over  several  scans.  This  can  be  represented  in  the  model 
by  the  following  procedure: 
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1 . Test  the  aircraft  turn  rate  for  a maneuver.  A suggested 

criterion  is  a turn  rate  magnitude  of  2 deg/sec  or  greater. 

2.  If  the  maneuver  test  is  positive,  generate  a random  number 
from  a uniform  distribution  between  0 and  1 . A lost  reply 

is  assumed  if  this  number  is  less  than  the  specified  probability 
for  signal  dropout.  The  probability  of  dropout  during  a turning 
maneuver  is  estimated  to  be  about  0.05  to  0.10. 

3.  Once  a particular  target  has  suffered  a dropout  in  step  2, 
increase  its  probability  of  dropout  during  the  remainder  of 
that  maneuver  to  about  0.8  to  0.9. 


The  output  of  the  ground  surveillance  model  is  shown  as  measured  range 
and  azimuth,  which  are  converted  to  rectangular  position  coordinates  in  the  radar 
tracking  model.  However,  the  conversion  could  occur  in  the  ground  surveillance 
model,  in  which  case  the  inputs  to  the  tracking  model  would  be  and  y^. 

The  rotating  characteristic  of  the  surveillance  antenna  has  not  been  included  in 
the  model,  i.e,  the  measured  positions  and  times  do  not  depend  upon  the  targets'  relative 
azimuths  from  the  antenna.  Instead,  the  model  provides  a snapshot  of  all  measured  air- 
craft positions  at  the  same  time  once  each  scan.  This  provides  a substantial  simplification 
in  the  model  implementation,  without  significantly  compromising  its  performance. 


Limits  (See  Fig . 22) 

Radar  is  subject  to  a maximum  elevation  anqle  y , which  restricts 
' max 

the  altitude  of  the  aircraft  in  the  vicinity  of  the  radar  for  a useful  return; 


h < h + 6,076. 1 "%/(x  - x )^  + (y  - y )^  tan  y (54) 

a r y a r v/a  1 r max 

The  norma  I elevation  angle  limit  is  30°  for  radar,  and  45°  for  ATCRBS  (Ref.  10). 

For  a usable  return,  the  aircraft  must  be  within  I ine-of-sight  of  the  radar 
antenna.  For  a circular  earth,  this  implies 


h > 

a 


0 


.788  yTxa  - *r)2  + (ya  - yr)2  - n/T 


(55) 
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Figure  22.  Ground  Survei  I lance  Model  Limits. 


The  slant  range  between  the  aircraft  and  ground  antenna  must  be  less  than 

R to  ensure  adequate  signal  strength  of  the  return.  R is  about  60  nm  for  ASR, 
max  3 max  ' 

and  around  200  nm  for  ARSR  and  ATCRBS.  The  radar  may  be  subject  to  a maximum 

altitude  restriction  h < h . For  ASR,  h is  specified  as  10,000  ft. 

max  max  r 

Implementation  of  the  model  could  be  slightly  simplified  by  not  including 
the  limits  explicitly.  In  this  case,  the  model  user  should  be  aware  of  the  restrictions 
and  make  sure  he  does  not  violate  them,  or  at  least  uses  them  to  interpret  his  results 
correctly. 


6.2  RADAR  TRACKING 

The  radar  tracker  provides  estimates  at  each  scan  of  the  target  position 
and  groundspeed  based  on  the  noisy  radar  data.  The  present  ATC  system  uses  an 
"a  - 3 tracker"  (Refs.  51  . 52)  which  provides  a weighted  first-order  correction 
to  the  position  and  groundspeed  each  scan.  More  sophisticated  systems  (Refs.  53  - 57) 
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have  been  investigated  and  could  provide  better  tracking  at  the  expense  of  more 
computation.  Figure  23  defines  important  geometry  and  nomenclature  for  the 
cr  - (3  tracking  model. 

The  measured  slant  range  and  azimuth  of  the  target  on  scan  i are  resolved 
into  approximate  rectangular  components: 

x.  = x + R sin  9 (56a) 

i r r.  r.  v ' 

i i 

y;  = yr  + Rr.  cos  er.  (56b) 

i i 


The  measured  position  coordinates  then  are  used  to  correct  the  predicted  position 
at  t.  which  was  calculated  on  the  previous  scan: 

x.  = x + o/(x.  - x ) 

i p.  i p/ 

ri  ri 

y.  = y + a{y.  - y ) 

'i  'p.  'i  'p. 

ri  ri 


(57a) 

(57b) 


The  parameter  a weights  the  contribution  of  the  measurement  to  the  updated  position 
estimate . 


The  estimated  groundspeed  components  from  the  previous  scan  are  updated, 
for  the  current  scan  using  x.  and  y.: 


a /% 


+ i(x. 


(58a) 


/S 


(58b) 


where,  in  the  first  scan  of  a new  target,  Xq  ~ 9 q = 0 • P is  a weighting  parameter 
which  determines  how  much  the  updated  groundspeed  changes.  The  estimated 
groundspeed  in  knots  is  computed  as 
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Figure  23.  Radar  Tracking  Geometry  and  Nomenclature. 
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(59) 


V = 3600 
9; 


* 


* 2 A 2 

+ y; 


The  current  estimated  position  and  groundspeed  components  are  then  used  to 
predict  the  aircraft  position  on  the  next  scan: 


P;+i 


^ i ^ t 

- x.  + x.T 

i i 


(60a) 


y 


Pi+1 


(60b) 


As  mentioned  in  the  Ground  Surveillance  Model  comments,  the  conversion 
from  polar  to  rectangular  coordinates  could  be  performed  in  either  model.  This  con- 
version has  been  included  in  the  tracking  model  for  two  reasons: 

a)  If  Mode  C (i.e.  automatic  altitude  reporting)  capability  is  available, 
the  tracking  computer  can  use  the  aircraft  altitude  to  correct  for 
slant  range  in  converting  to  rectangular  coordinates.  At  high 
elevation  angles,  the  slant  range  can  introduce  position  errors 

of  more  than  40% . 

b)  In  the  ARTS  111  tracker,  radar  range  and  azimuth  are  used  directly 
in  correlating  target  returns  with  aircraft  tracks.  Thus,  if  a more 
detailed  model  of  the  ARTS  III  tracker  is  developed  later,  it  will 
be  compatible  with  the  present  models. 


A major  element  of  the  ARTS  III  tracker  (Ref.  58)  is  the  correlation  of 
beacon  replies  with  active  tracks.  Because  the  beacon  reply  is  frequently  lost, 
garbled,  reflected  or  duplicated,  ARTS  III  includes  an  elaborate  and  complicated 
logic  for  correlation.  The  confidence  level  of  the  track's  position  is  a function 
of  the  number  of  successful  correlations.  Whenever  a track  is  not  successfully 
correlated,  its  confidence  level  is  decreased,  and  the  estimated  position  is  equated 
to  its  predicted  position.  On  a successful  correlation,  the  predicted  position  is 
updated  by  Equation  (57)  and  the  confidence  level  increased. 

The  magnitudes  of  the  update  parameters  (»,p)  are  functions  of  the  individual 
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track's  confidence  level:  the  higher  the  confidence  level,  the  more  weight  is 
placed  on  its  predicted  position  and  the  less  on  the  new  measurement  (see  Table  3). 

Table  3.  Variation  of  ARTS  III  Tracking-Parameters  with 
Successful  Correlation. 
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At  the  beginning  of  a track  the  position  measurement  is  used  entirely  and  ground- 
speed  is  not  calculated.  After  the  track  has  successfully  correlated  for  11  or  more 
scans,  the  weighting  parameters  reach  a steady-state  condition.  On  an  unsuccessful 
correlation,  the  sizes  of  the  weighting  parameters  are  increased  slightly  (about 
1/2  step).  The  track  is  placed  in  a permanent  coast  mode  and  the  controller  alerted 
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after  three  successive  unsuccessful  correlations. 

Two  versions  of  the  tracker  could  be  implemented  in  the  TSC  simulation 
without  requiring  the  very  complex  correlation  logic  in  ARTS.  The  simpler  tracker 
would  use  constant  values  of  a and  |3;  the  steady-state  values  in  Table  3 are  recom- 
mended . A more  accurate  model  would  include  the  variation  of  a and  (3  shown  in 
Table  3.  In  the  event  of  lost  replies  the  track  would  be  coasted,  as  discussed  above, 
but  cr  and  (3  would  remain  at  their  previous  values  since  the  actual  increases  are 
very  smal  I . 

The  correlation  and  updating  calculations  of  the  tracker  require  a finite  period 
of  time  to  perform.  In  ARTS  III,  the  position  and  groundspeed  estimates  are  delayed 
by  about  one  second  after  the  surveillance  measurements.  This  delay  is  probably  not 
significant  enough  to  warrant  inclusion  in  the  model. 

As  in  the  surveillance  model,  the  effects  of  the  relative  azimuths  of  the  targets 
are  not  included  in  the  tracking  model.  The  tracks  are  all  updated  at  the  same  time 
each  scan.  This  eliminates  such  complications  as  modifying  the  predicted  position  to 
account  for  a change  in  measurement  time  due  to  a change  in  azimuth. 
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SECTION  7 


DISCUSSION 

The  models  presented  in  the  previous  sections  have  been  developed  to  provide 
more  accuracy  and  realism  for  simulations  requiring  the  high  frequency  response  of 
advanced  air  traffic  control  systems.  Addition  of  these  models  to  the  existing  program 
will  add  greater  realism  to  the  simulation,  but  will  also  require  more  computer  time. 

In  general,  the  more  sophisticated  models  should  only  be  used  when  the  simpler  models 
are  inadequate.  The  following  discussion  compares  the  results  one  would  expect  using 
the  simple  as  opposed  to  the  complex  model  and  indicates  when  it  would  be  desirable 
to  use  the  more  realistic  model. 

The  ROSS  simulation  modeled  the  wind  by  layers  of  constant  direction  and 
velocity.  Early  In  the  present  study  it  was  recommended  that  a more  reasonable  model 
would  provide  a linear  variation  of  both  wind  speed  and  direction  with  altitude.  This 
suggestion  has  since  been  implemented  in  the  simulation  program.  It  was  concluded 
the  effects  of  vertical  winds  and  higher  frequency  wind  gusts  would  be  much  smaller 
than  the  other  errors  being  modeled,  and  therefore  unnecessary  to  include. 

In  the  present  program  the  airspeed  does  not  change  immediately  with 
wind  changes  but  the  groundspeed  does  suddenly  jump  to  its  steady-state  value. 

The  real  aircraft  responds  to  sudden  wind  changes  by  having  its  airspeed  suddenly 
change  to  the  new  steady-state  value  with  the  groundspeed  eventually  catching 
up,  an  effect  which  is  included  in  the  more  sophisticated  model.  Normally  the 
difference  between  the  two  models  is  insignificant  because  the  pilot  is  using 
thrust  to  hold  airspeed  constant  in  a fairly  tight  control  loop.  The  more  sophisticated 
model  would  be  desirable  when  one  is  using  groundspeed  instead  of  airspeed  commands 
such  as  in  4-D  guidance,  or  when  there  is  concern  about  the  effects  of  gusts  or 
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severe  wind  shears. 

Altitude  in  the  present  program  is  changed  at  a constant  rate  until  it 
reaches  its  new  value,  where  it  stops  abruptly.  The  new  model  uses  two  rates  and 
makes  the  transition  more  realistic  by  including  a response  lag  and  by  limiting  the 
normal  acceleration.  The  more  sophisticated  model  would  be  desirable  where 
altitude  separations  are  critical,  such  as  the  evaluation  of  CAS  or  IP.C.  Figures 
24  - 27  compare  the  altitude  response  models  for  four  aircraft  ranging  from  single- 
engine propeller  to  standard  commercial  jet.  The  ASI  model  results  were  obtained 
by  digital  simulation  of  the  model  described  in  Section  3.1,  neglecting  any  air 
data  system  errors.  In  the  revised  (ASI)  model,  the  rate  of  climb  or  descent  decreases 
to  500  fpm  when  the  altitude  difference  is  less  than  1,000  ft.  Because  of  the  altitude 
tolerance  in  the  new  model,  the  aircraft  begins  to  level  off  when  it  is  within  50  ft 
of  the  desired  altitude.  The  simulation  was  stopped  when  the  altitude  difference 
was  within  50  ft  and  the  altitude  rate  was  less  than  50  fpm. 

Airspeed  changes  in  the  present  system  are  also  at  a constant  rate  until 
reaching  the  desired  value.  The  new  model  introduces  a time  lag  into  the  application 
of  acceleration,  and  also  limits  acceleration  during  climb  and  deceleration  during 
descent  in  accordance  with  physical  reality.  The  more  sophisticated  model  would 
be  desirable  particularly  where  velocity  control  is  being  used  for  separation,  spacing 
or  merging.  Figures  28  - 31  compare  the  responses  of  the  two  models  to  commanded 
airspeed  changes,  simultaneously  with  the  altitude  changes  shown  in  Figures  24  - 27. 
Again  the  ASI  results  were  obtained  by  simulation,  neglecting  the  air  data  system 
errors  and  assuming  no  winds.  The  simulation  was  terminated  when  the  airspeed 
error  was  less  than  2.5  kts  and  the  acceleration  (or  deceleration)  was  within  10% 
of  ACR  (or  DCR),  and  the  altitude  and  altitude  rate  criteria  above  were  satisfied 
simultaneously.  The  airspeed  response  of  the  present  model  is  unaffected  by  the 
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ALTITUDE  (h),  FT 


Figure  24.  Comparison  of  ROSS  and  ASI  Altitude  Response  Models 
for  a Climb  and  a Descent  (Single-Engine  Propeller  Air- 
craft) . 
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ALTITUDE  (h),  FT 


15,000 


Figure  25.  Comparison  of  ROSS  and  ASI  Altitude  Response  Models 

for  a Climb  and  a Descent  (Medium  Turbo-Prop  Aircraft). 
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ALTITUDE  (h), 


25,000 


Figure  26.  Comparison  of  ROSS  and  ASI  Altitude  Response  Models 
for  a Climb  and  a Descent  (Executive  Jet  Aircraft). 
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ALTITUDE  (h),  FT 


Figure  27.  Comparison  of  ROSS  and  ASI  Altitude  Response  Models 
for  a Climb  and  a Descent  (Standard  Commercial  Jet 
Aircraft) . 
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AIRSPEED  (Va),  KT 


Figure  28.  Comparison  of  ROSS  and  ASI  Velocity  Response  Models 
for  Climbing  Speed-Up  and  Descending  Slow-Down 
(Single  Engine  Propeller  Aircraft). 
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AIRSPEED  (Va),  KT 


Figure  29.  Comparison  of  ROSS  and  ASI  Velocity  Response  Models 
for  Climbing  Speed-Up  and  Descending  Slow-Down 
(Medium  Turbo-Prop  Aircraft). 
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Figure  30.  Comparison  of  ROSS  and  ASI  Velocity  Response  Models 
for  Climbing  Speed-Up  and  Descending  Slow-Down 
(Executive  Jet  Aircraft). 
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AIRSPEED  (Va),  KT 


Figure  31 . Comparison  of  ROSS  and  ASI  Velocity  Response  Models 
for  Climbing  Speed-Up  and  Descending  Slow-Down 
(Standard  Commercial  Jet  Aircraft) . 
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altitude  response.  However,  in  the  revised  model  the  altitude  response  limits 
the  acceleration  or  deceleration  of  the  aircraft.  In  the  slowing-down,  descending 
cases,  the  limit  is  not  as  effective  since  the  normal  rate  of  descent  is  somewhat 
less  than  the  maximum  (see  Table  1). 

Heading  changes  also  occur  at  constant  rate  in  the  present  program. 

The  more  sophisticated  model  allows  a time  lag  for  establishing  the  turn  plus  a 
simulation  of  the  pilot's  turn  logic  which  may  cause  slight  over  or  under  shoot 
of  the  desired  heading.  The  more  sophisticated  model  would  be  desirable  where 
turn  maneuvers  are  used  for  CAS  or  IPC,  and  also  where  path  stretching  is  used 
for  separation,  spacing  or  merging.  The  two  model's  responses  are  compared  in 
Figures  32  - 34  for  heading  changes  of  10°,  30°  and  90°  respectively.  The  simulation 
results  in  Figures  32  - 34  are  for  constant  altitude  and  airspeed,  and  neglect  any 
air  data  system  errors.  The  simulation  was  terminated  when  the  heading  was 
within  0.5°  of  the  commanded  value  and  the  turn  rate  was  less  than  0.1  °/sec. 

Aircraft  guidance  along  a desired  track  in  the  present  simulation  is 
assumed  to  be  perfect  after  the  aircraft  actual  position  is  within  1 .5°  of  the  desired 
bearing.  The  new  model  for  guidance  along  a desired  track  is  much  more  sophisticated 
It  uses  a steering  algorithm  which  depends  on  the  distance  off  track  ana  the  heading 
relative  to  the  desired  track  angle.  In  order  to  acquire  the  track  the  aircraft  is 
first  placed  on  an  intercept  heading  of  approximately  90°  to  the  desired  course. 

Then  when  sufficiently  close,  the  algorithm  is  changed  to  that  of  a linear,  second- 
order  regulator  which  causes  the  aircraft  to  acquire  the  track  with  a slight  overshoot. 
Steady-state  offset  can  occur  because  of  a crosswind.  The  logic  is  similar  to  what 
is  used  in  coupled  automatic  pilots.  The  more  sophisticated  model  would  be  preferred 
whenever  the  acquisition  of  an  assigned  track  is  left  to  the  pilot  or  whenever  the 
tracking  dynamics  are  of  importance  such  as  with  closely  spaced  parallel  track 
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Figure  32.  Comparison  of  ROSS  and  ASI  Heading  Response  Models 
for  Various  Aircraft  Categories  (10°  Heading  Change). 
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Figure  33.  Comparison  of  ROSS  and  ASI  Heading  Response  Models 
for  Various  Aircraft'  Categories  (30°  Heading  Change). 
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Comparison  of  ROSS  and  ASI  Heading  Response  Models 
for  Various  Aircraft  Categories  (90°  Heading  Change). 
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airways  or  ILS  installation. 

The  navigation  error  in  the  present  program  consists  of  correlated  noise 
added  to  the  true  cartesian  position  coordinates.  Both  components  of  position  noise 
use  the  same  standard  derivation  and  correlation  time.  Navigation  errors  in  the 
new  models  are  more  representative  of  specific  navigation  systems  which  might 
be  in  use.  Various  error  sources  can  be  identified  to  ascertain  relative  importance 
of  errors.  Geometric,  computation  and  propagation  errors  can  be  changed  independ- 
ently. The  more  sophisticated  model  would  be  desirable  when  the  objective  of 
the  experiment  is  to  ascertain  the  relative  importance  of  different  types  of  error 
sources,  or  to  evaluate  competing  navigation  systems  for  particular  applications. 

The  present  surveillance  error  is  computed  by  random  selection  of  separate 
range  and  bearing  bias  errors  defined  by  mean  and  standard  deviation  input  values.  The 
new  model  changes  the  range  and  bearing  errors  on  each  discrete  scan  and  processes  the 
radar  data  using  the  same  basic  tracking  algorithms  used  in  the  third-generation  system. 
The  more  sophisticated  model  should  be  used  whenever  the  radar  data  is  to  be  processed 
for  velocity  determination,  conflict  prediction,  radar  mosaics,  etc. 

It  should  be  clear  that  the  more  sophisticated  models  can  be  used  selectively, 
one  at  a time,  depending  upon  the  specific  need.  They  might  be  applied  only  to  one  or 
two  aircraft  in  a simulation  or  they  can  be  invoked  only  when  two  aircraft  are  in  a con- 
flict situation.  Used  in  this  manner,  the  program  is  longer,  but  the  increase  in  running 
time  can  be  kept  minimal . Experience  with  the  program  should  show  when  the  more 
sophisticated  models  are  justified  and  when  their  increased  realism  is  superfluous. 

In  summary,  Table  4 presents  a comparison  of  the  present  and  revised 
models  for  the  more  Important  ATC  system  functional  error  sources.  As  indicated, 
the  effects  of  certain  of  the  lesser  important  or  indescribable  error  sources  have  not 
yet  been  included  in  either  model. 
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Table  4.  Comparison  of  ATC  System  Functional  Error  Models. 


FUNCTIONAL  STANDARD  REVISED 

ERROR  SOURCE  ROSS  MODEL  ASI  MODEL 


AIR  DATA  SYSTEM 

Altimeter 

★ 

random  bias 

bias  a function  of  h 

Rate  of  cl imb 

* 

random  bias 

bias  proportional  to  h 

Heading 

★ 

random  bias 

random  bias 

Airspeed/Mach  No. 

★ 

random  bias 

bias  proportional  to  V 

NAVIGATION  SYSTEM 
VOR 

DME 

correlated  noise 
in  cartesian 
coordinates  - 
no  directivity 

geometry,  propagation,  site 
& receiver  errors  modeled 

geometry,  propagation,  site 
& receiver  errors  modeled 

RNAV 

can  be  approximated 
as  perceived  x and  y 
location  errors  relative 
to  a described  location 

computed  from  VOR/DME 
errors 

ILS 

not  modeled 

not  modeled 

AIRCRAFT  DYNAMICS 

Acceleration 

★ 

no  limits 
instantaneous 

limited  by  h 
first  order  lag 

Climb  rate 

* 

single  rate 
instantaneous 

variable  rate 
first  order  lag 

The  description  given  applies  to  "standard"  use  of  the  simulator.  At  his  option, 
the  user  can  prepare  route  procedures  that  achieve  other  conditions.  For  example 
the  bias  of  a value  can  be  made  proportional  to  the.  value,  one  rate  can  be  made 
a function  of  another,  etc. 
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Table  4.  (Cont.) 


FUNCTIONAL  STANDARD  REVISED 

ERROR  SOURCE  ROSS  MODEL  ASI  MODEL 


AIRCRAFT  DYNAMICS  (Cont.) 

Turn  rate 

* 

single  rate 

variable  rate 

instantaneous 

lag  for  roll-in 

Tracking 

acquisition  threshold 

complete  acquisition  dynamics 

no  tracking  error 

2nd  order  linear  tracker 

Pilot  reaction  time 

not  modeled 

included  in  lags 

Blunders 

not  modeled 

not  modeled 

SURVEILLANCE  SYSTEM 

Range  error 

fixed  bias 

random  bias  each  scan 

Azimuth  error 

fixed  bias 

random  bias  each  scan 

Geometry  error 

not  included 

included 

Tracking  algorithm 

not  included 

simulated  actual  algorithm 

WIND 

Variation  with  altitude 

steady-state  only 

proper  aircraft  response 

Variation  with  location 

none 

proper  aircraft  response 

Prediction  error 

not  modeled 

not  modeled 

The  description  given  applies  to  "standard"  use  of  the  simulator.  At  his  option, 
the  user  can  prepare  route  procedures  that  achieve  other  conditions.  For  example 
the  bias  of  a value  can  be  made  proportional  to  the  value,  one  rate  can  be  made 
a function  of  another,  etc. 
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Table  4.  (Cont.) 


FUNCTIONAL  STANDARD  REVISED 

ERROR  SOURCE  ROSS  MODEL  ASI  MODEL 


COMMUNICATIONS 

Delay 

Garble 


not  modeled 
not  modeled 


not  modeled 
not  modeled 


CONTROLLER 
Reaction  time 
Decision  delay 
Blunders 


not  modeled 
not  modeled 
not  modeled 


not  modeled 
not  modeled 
not  modeled 
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SECTION  8 


CONCLUSIONS  AND  RECOMMENDATIONS 

This  section  presents  a summary  of  the  important  conclusions  of  the  error 
analysis  and  modeling  study,  along  with  some  suggested  applications  for  the  revised 
error  models  and  possible  additional  improvements  for  the  simulation. 

8.1  SUMMARY  OF  CONCLUSIONS 

The  important  functional  error  sources  for  the  present  and  near-future 
ATC  system  were  identified  and  analyzed  for  their  relative  effects  on  system 
performance.  Error  models  were  developed  for  the  dominant  functional  error  sources 
having  response,  times  of  the  order  of  three  seconds  or  greater.  The  new  models 
thus  provide  approximately  an  order  of  magnitude  improvement  for  the  ATC  simulation 
over  the  existing  models  which  generally  neglect  responses  of  less  than  about  30 
seconds. 

The  new  models  provide  the  simulation  user  with  the  capability  of  experi- 
menting with  a number  of  future  ATC  system  concepts  which  could  not  be  properly 
evaluated  with  the  present  programs.  Those  ATC  concepts  which  depend  on  higher 
frequency  and  more  accurate  responses  can  be  effectively  tested  with  confidence 
that  the  major  error  sources  are  included.  Examples  of  the  experiments  which  can 
be  investigated  with  the  new  models  include  collision  avoidance  systems,  intermittent 
positive  control,  4-D  navigation,  and  metering  and  spacing  concepts. 

The  target  dynamics  models  have  been  developed  to  duplicate  as  nearly 
as  practicable  the  normal  response  that  an  actual  pilot  would  produce.  This 
required  tradeoffs  between  computational  complexity  and  realism.  The  models  were 
constructed  to  be  as  independent  of  specific  aircraft  performance  data  as  feasible. 
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As  a result,  the  representation  of  a particular  aircraft  requires  the  specification 
of  sixteen  or  fewer  parameters  and  time  constants.  The  altitude  and  airspeed/ 
heading  control  loops  are  uncoupled,  except  that  the  acceleration-deceleration 
capability  is  realistically  limited  during  a simultaneous  climb-descent. 

The  altitude,  airspeed  and  heading  control  models  contain  nonlinear 
logic  elements  which  required  some  verification  of  their  intended  operation. 

Simple  check  programs  were  developed,  neglecting  the  air  data  system  and  wind 
effects,  to  evaluate  the  performance  of  these  models.  The  results  of  the  check 
programs  conform  to  the  desired  response. 

Analysis  of  the  wind  effects  on  system  performance  indicated  that  gusts 
will  not  be  significant  for  the  advanced  ATC  concepts  mentioned  previously. 
Consequently,  gusts  were  neglected  in  the  revised  wind  model.  However,  unknown 
steady-state  winds  and  strong  wind  shears  do  appear  to  present  problems,  particularly 
for  4-D  guidance  and  metering  and  spacing  concepts. 

The  or  - (3  tracker  algorithm  in  the  surveillance  model  introduces  a significant 
delay  in  sensing  groundspeed  and  can  result  in  substantial  errors  in  the  estimated 
target  position.  The  tracker  was  originally  designed  for  the  purpose  of  keeping 
the  alphanumerics  on  the  radar  plan  position  indicator  associated  with  the  correct 
target;  it  was  not  intended  to  estimate  groundspeed  of  the  target.  The  a - 3 tracker 
is  based  on  nonaccelerating  flight,  and  consequently  cannot  accurately  respond 
to  rapid  target  maneuvers. 

8.2  RECOMMENDED  APPLICATIONS 

Following  are  some  recommended  applications  for  the  ATC  simulation 
model.  In  most  cases  the  experiments  can  be  conducted  with  existing  software 
and  the  newly-developed  models  presented  in  this  report. 
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1 . Terminal  Area  Simulation  - Error  Assessment. 

The  purpose  of  this  experiment  is  to  examine  the  effects  of  wind,  air 
data  system  errors,  navigation  errors  and  surveillance  errors  on  safety,  capacity, 
and  control  effort  in  the  terminal  area.  The  ATC  simulation  is  exercised  using 
a predetermined  control  process,  while  the  relative  contributions  of  the  individual 
error  sources  are  varied  systematically.  Safety  is  measured  by  the  number  of  times 
a given  separation  criterion  is  violated.  Capacity  is  defined  by  the  number  of 
aircraft  handled  per  unit  time.  Control  effort  can  be  estimated  by  the  number 
of  control  commands  which  are  necessary. 

2.  Enroute  Area  Simulation  - Error  Assessment. 

The  purpose  of  this  experiment  is  to  examine  the  effects  of  wind,  air 
data  system  errors,  navigation  errors  and  surveillance  errors  on  the  accuracy  of 
metering  aircraft  to  the  terminal  area.  The  ATC  simulation  is  programmed  to  feed 
aircraft  to  the  terminal  area  landing  process  while  each  of  the  error  sources  is 
varied  systematically.  Metering  accuracy  is  measured  by  the  timing  error  at  the 
terminal  area  input  fix. 

3.  4-D  Guidance  Simulation  - Error  Assessment . 

The  purpose  of  this  experiment  is  to  examine  the  effects  of  wind,  air 
data  system  errors,  navigation  errors  and  surveillance  errors  on  the  ability  of  a 
target  to  follow  a 4-D  path.  A typical  arrival  path  must  be  established  with  arrival 
fixes  designated  along  it.  Then,  a control  logic  for  generating  discrete  commands 
to  the  aircraft  must  be  developed.  These  commands  should  include  both  speed 
control  and  path  stretching  in  order  to  correct  for  disturbances  which  change  the 
arrival  time.  The  accuracy  of  time  of  arrival  will  be  measured  in  the  face  of  wind, 
air  data,  navigation  and  surveillance  system  errors.  Knowledge  of  an  unsteady 
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wind  is  expected  to  be  an  important  factor  and  the  simulation  should  establish 
how  accurately  the  wind  should  be  known  and  to  what  granularity  as  a function 
of  space  and  time  . 

4.  Inter-Aircraft  Spacing  Assessment. 

A spacing  and  sequencing  process  which  uses  actual  instead  of  planned 
information  on  the  positions  of  targets  may  have  efficiency  and  capacity  advantages 
over  fixed  slot,  preplanned  processes.  Previous  applications  have  modeled  the  single 
aircraft  when  it  is  supplied  a flight  plan  and  position  relative  to  navigational 
coordinates,  but  the  important  questions  for  the  total  system  concern  the  inter- 
aircra-ft  spacings.  The  strategy  taken  by  both  the  controller  and  the  individual 
aircraft  commander  must  be  modified  in  that  each  may  now  be  influenced  by  the 
presence  of  additional  aircraft.  This  requires  development  of  a new  sequencing 
and  spacing  algorithm,  including  the  insertion  of  take-offs  in  gaps  which  open  up 
in  the  ianding  flow. 

Interactive  simulation  using  human  pilots  and  controllers  is  recommended 
first  inasmuch  as  the  existing  system  does  not  provide  traffic  information  in  the 
cockpit  concerning  nearby  aircraft.  Such  a simulation  does  not  need  to  be  overly 
complex;  each  pilot  needs  only  three  controls  to  set  velocity  rate,  heading  rate 
and  altitude  rate.  Aircraft  dynamics  and  constraints  described  by  the  AS1  models 
can  be  used  to  simulate  various  aircraft.  The  positions  of  all  the  simulated  aircraft 
can  be  presented  on  a common  display.  After  pilot  strategy  has  been  determined 
by  observation  it  can  be  modeled  for  analytic  simulation.  The  controller  may  be 
simulated  by  either  a human  or  an  analytic  algorithm.  Distributed  tactical  systems 
may  be  simulated  by  letting  the  human  pilots  operate  without  the  conttoller  in 
the  loop.  IRC  systems  may  be  simulated  by  modification  of  the  strategy  used 
by  the  controller  and  eliminating  the  pilot's  access  to  the  display.  The  first  task 
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is  to  develop  the  interactive  simulation.  The  second  task  is  to  observe  the  inter- 
actions and  model  the  strategies  used  by  both  pilots  and  controllers  under  the  various 
control  concepts.  The  next  task  is  to  exercise  the  analytic  models  and  compare 
the  results  to  those  obtained  with  human  pilots  and  controllers.  Finally,  the 
analytic  models  can  be  exercised  to  evaluate  the  various  control  concepts.  Capacity 
and  risk  evaluations  can  be  determined  from  the  inter-aircraft  spacings  which  result. 

5.  Separation  Assurance  Concepts  Evaluation. 

The  ATC  simulation  can  be  applied  to  the  evaluation  of  proposed  CAS 
and  IPC  concepts.  At  present  the  traffic  situation  can  be  programmed  so  that 
targets  fly  prescribed  tracks  and  respond  to  input  commands.  A number  of  new 
concepts  and  codings  will  have  to  be  developed.  First,  a hazard  criterion  must 
be  defined  in  terms  of  relative  position,  relative  speed  and  relative  acceleration. 
Algorithms  must  be  established  for  checking  collision  pairs  for  violation  of  the 
hazard  criterion,  and  for  resolving  the  hazard  and  issuing  appropriate  commands 
to  the  aircraft  involved.  Also,  the  aircraft  has  to  be  returned  to  normal  flight 
following  the  hazard  resolution.  Finally,  coding  must  be  provided  for  keeping 
statistics  on  true,  false  and  missed  alarms,  duration  of  hazards,  distance  of  closest 
approach,  and  other  evaluation  parameters. 
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